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INTRODUCTION 


This document describes the Langley Research Center 16-Foot Transonic Tunnel standard 
set of equations. The engineering units necessary for these equations are computed on site 
from the raw data millivolts or counts. These quantities with additional constants are used 
as input to the program for computing the forces and moments and the various 
coefficients. 

This document is intended to be a companion document to NASA Technical Memorandum 
83186, A User's Guide to the Langley 16-Foot Transonic Tunnel, August 1981. 

The equations are grouped into modules, so that only the required modules need be used. 
The modules are as follows: 

A. Wind Tunnel Parameters 

B. Jet Exhaust Measurements' 

C. Skin Friction Drag 

D. Balance Loads and Model Attitudes 

E. Internal Drag (or Exit-Flow Distributions) 

F. Pressure Coefficients and Integrated Forces 

G. Thrust Removal Options 

II. Turboprop Options 

Individual customizing of these equations for a specific job application is permitted through 
the use of code constants. These equations do not cover all possible jobs; however, they 
are coded so that modifications of selected equations may be easily carried out. 



The format of this document is arranged so that the module designations correspond to the 
Appendix designations in which the respective calculations equations are given. 

WIND TUNNEL PARAMETERS 

The wind tunnel parameters are computed from the required static and total pressure 
measurements. The Reynolds number, dynamic pressure and tunnel total temperatures are 
computed. When the tunnel Mach number is computed, a lookup table from an earlier wind 
tunnel calibration is used to correct the ratio of static pressure to total pressure used in 
the Mach number calculation. These wind tunnel parameters are stored for use by other 
modules. Refer to Appendix A for calculations. 

JET EXHAUST MEASUREMENTS 

Jet exhaust information is calculated for the primary, secondary and tertiary flow 
conditions. 


The primary flow conditions for each engine, up to a maximum of four, are calculated. The 
various parameters that are computed are mass flow and ideal thrust for each engine. The 
average nozzle pressure ratio and average total temperature over all the engines is 
obtained. The total mass flow is derived from chamber and/or venturi measurements. 
Discharge coefficients for the total system are computed as well as the ideal thrust. 

For the secondary and tertiary flows, the mass flows and other parameters are computed. 
Refer to Appendix B for calculations. 

SKIN FRICTION DRAG 

The skin friction drag for the model is computed in addition to any empennage skin friction 
drag. Refer to Appendix C for calculations. Information from the wind tunnel parameters 
is used. Drag from the various components as well as total drag is computed. 
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BALANCE LOAD AND MODEL ATTITUDES 


The balance computations for the force and moment coefficients for up to three balances 
may be computed from this module. Allowances for the method of attaching the balances 
are made. The measured forces and moments are corrected for balance interactions. Then 
an allowance is made for high order interactions and momentum tares. The forces and 
moments are rotated to the desired axis and the final correct coefficients are computed as 
well as the angle of attack and sideslip angles. Refer to Appendix D for calculations. 

INTERNAL DRAG (or EXIT-FLOW DISTRIBUTIONS) 


The internal drag and various forces on the engines are computed using the equations given 
in Appendix E. The result of these computations are used in the balance computations of 
module D to correct the force measured by the balances. 

PRESSURE COEFFICIENTS AND INTEGRATED FORCES 

Pressure coefficients are computed by using the equations given in Appendix F. Various 
integrated forces due to the pressures are calculated including hinge moment coefficients. 

THRUST REMOVAL 


Various thrust removal coefficients may be computed according to specified flags which 
specify the model setup. Various configurations are permitted which may include two 
balances. Reference Appendix G for calculations. 

TURBOPROP OPTIONS 

The drag and thrust coefficients due to the propeller and jet engine are computed as well 
as the combined totals. Horsepower and efficiency of the engines are derived with other 
quantities. Reference Appendix H for calculations. 
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APPENDIX A 
Tunnel Parameters 


Nomenclatures A-l 

Required Constants A -3 

Atmospheric Pressure A-3 

Mach Number A-3 

Tunnel Static Pressure A-4 

Tunnel Total Pressure A-5 

Tunnel Dynamic Pressure A-5 

Dew Point A -6 

Tunnel Total Temperature A-6 

Reynolds Number A-6 
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MODULE A TUNNEL PARAMETERS 


A-l 


SYMBOL 

MACH 

MCODE 


PO 

PO/PTO 

PTANKG 

PTANKH 

PTG 

PTH 

PTKSON 

PTO 

PTSON 

QO 

REFL 

RN 

RN/FT 

RT(J) 

T(J) 

TTO 


NOMENCLATURE 
Free stream Mach number. 

Mach number calculation code. 

=1, PTANKG and PTH are needed. 

=2, PTANKH and PTH are needed. 

=3, PTANKG and PTG are needed. 

=4, PTANKH and PTG are needed. 

=5, PTKSON and PTSON are needed. 

Tunnel static pressure, lbs/sq. in. 

Ratio of tunnel static pressure to total pressure. 

Tunnel tank pressure measured by gage, lbs/sq. in. 

Tunnel tank pressure measured by Ruska, lbs/sq. in. 

Tunnel total pressure measured by gage, lbs/sq. in. 

Tunnel total pressure measured by Ruska, lbs/sq. in. 

Tunnel tank pressure measured by Digiquartz, lbs/sq. in. 
Tunnel total pressure, lbs/sq. in. 

Tunnel total pressure measured by sonar manometer, lbs/sq. 
in. 

Dynamic pressure, lbs/sq. in. 

Reference length, feet. 

Reynolds number based on reference length. 

Reynolds number per foot. 

Tunnel total temperature measurements, °F, 
where J = probe number. 

Constants required from project engineer (0.0 or 1.0) 
where J = probe number. 

Tunnel total temperature, °F. 
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APPENDIX A 
Module A 

Tunnel Parameters 


Required Constants 

1. MCODE (default value = 2) must be provided if values other than PTKSON 
and PTSON are used to compute Mach number. 

2. The constants used in determining tunnel total temperature are T2, T3, T4 
and T5 which must equal 0.0 or 1.0. 

One-tunnel temperature measurement 

T2 = 1.0, T3 = T4 = T5 = 0.0 (Eq- A-l) 

Two-tunnel temperature measurements 

T2 = T3 = 1.0, T4 = T5 = 0.0 (Eq. A-2) 

Note that the numbers 2 through 5 correspond to resistance thermometer 
numbers normally used. 

3. A reference model length, REFL, must be given in units of feet to 
compute model Reynolds number. 

Atmospheric Pressure 

Atmospheric pressure calculation may be handled in the standard program for 
quantities. Its inclusion (if required) and method of obtaining (dialed-in optional 
digital channel or measured by gauge in analog channel) is left optional to the 
project engineer. However, measuring atmospheric pressure with a gauge is 
recommended rather than entering this pressure reading into an analog channel 
since it is possible for significant variations to occur during the course of a 
tunnel run. 

Mach Number 

1. MCODE indicates which measurements are to be used for Mach number 
calculation (see nomenclature on page A-l). The default value of MCODE 
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A-4 

is 2. Multiple options are provided to allow for the possibility of 
instrument failure during a test. If the digital MCODE input is 1 to 5, 
then digital value overrides the C-card value. If the digital value is zero, 
then the "C” value overrides. The reference pressures may also change. 

If MCODE = 1 

PO/PTO = (PTANKG/PTH)K + I (Eq. A-3) 

If MCODE = 2 

PO/PTO = (PTANKH/PTH)K + I (Eq. A-4) 

If MCODE = 3 

PO/PTO = (PTANKG/PTG)K + I (Eq. A-5) 

If MCODE = 4 

PO/PTO = (PTANKH/PTG)K + I (Eq. A-6) 

If MCODE = 5 

PO/PTO = (PTKSON/PTSON)K + I (Eq. A-7) 

where K and I are from 1965 16-ft TT calibration 

MACH =75 ((PO/PTO)" 2/7 - 1) (Eq. A-8) 

D. Tunnel Static Pressure 

PO calculation automatically depends on MCODE. No input is required from 
the project engineer. The normal procedure (internal constant MCODE = 5) 
uses PTSON for computation. 
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A -5 

If MCODE < 2 

PO = (PO/PTO)PTH (Eq. A-9) 

If MCODE = 3 or 4 

PO = (PO/PTO)PTG (Eq. A-10) 

If MCODE = 5 

PO = (PO/PTO)PTSON (Eq- A-ll) 

E. Tunnel Total Pressure 

PTO to calculation automatically depends on MCODE. No input is required 
from the project engineer. The normal procedure (MCODE = 5) uses PTSON. 

If MCODE < 2 

PTO = PTH (Eq* A -12) 

If MCODE = 3 or 4 

PTO = PTG (Eq* A -13) 

If MCODE = 5 

PTO = PTSON (Eq- A-14) 


F. Tunnel Dynamic Pressure 

Tunnel dynamic pressure is computed as follows: 


If MACH < .1 


QO = PO 


(Eq. A-15) 


If MACH — .1 


QO = 0.7 * PO * MACH 2 


(Eq. A-16) 


Dew Point 

Dew point calculation may be handled in the standard program for quantities. 

Its inclusion, channel location, and name are left optional to the project 

engineer; however, TDP is suggested as a name. 

Tunnel Total Temperature 

1. Provision is made for four individual tunnel total temperature 
measurements. They may be either thermocouples or resistance 
thermometers; however, the appropriate equation must be specified for 
the standard program for quantities. Note that resistance thermometer 
one (1) (strut head) should not be used. If resistance therometers are 
used, their calibrations are included internal to the program. 

2. The constants required from the project engineer are T2, T3, T4, and T5 
(0.0 or 1.0). 


(RT2 * T2) + (RT3 * T3) + (RT4 * T4) + (RT5 * T5) 
TTO - T2 + T3 + T4 + T5 


(Eq. A -17) 


Reynolds Number 

1. The constant required from the project engineer is REFL. 
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RN = RN/FT * REFL 


(Eq. A-18) 
(Eq. A-19) 



APPENDIX B 


Jet Exhaust Measurements 


Nomenclatures 

B-l 

Required Constants 

B-9 

Test for Exhaust Model 

B-10 

Compute Common Constants 

B-10 

Individual Engine Measurements 

B-ll 

Total Exhaust System Properties 

B-l 6 

Secondary Flow Measurements 

B-21 

Tertiary Flow Measurements 

B-23 


MODULE B JET EXHAUST MEASUREMENTS 


B-l 


SYMBOL 

AENG(I) 


AREF 

AT(I) 

AVRI(L) 


C* 

CDSKL) 

CFI 

CFICHR 

FI 

FICHR 

FIENG(I) 


FM1 

FMS 

FVRI(I) 

GAMJ 

ICH(I) 


INTFM1 

INTFMS 

KAE(I) 


NOMENCLATURE 

Flow area to be used for determining each engine mass-flow 
rate from plenum chamber measurements, where I = engine 
number. This area is generally based on the area of the 
plenum orifice nozzles (AENG(I) = (orifice area)/2 for twin 
engines), sq. in. 

Model reference area used for coefficients, sq. in. 

Throat area of each engine, where I = engine number, sq. in. 
Area of throat of in-line (not MCV) venturi, where L = 
venturi number, sq. in. 

Critical area, sq. in. 

Discharge coefficient, where L = venturi number. 

Ideal thrust coefficient based on measured mass-flow rate. 
Ideal thrust coefficient based on mass-flow rate obtained 
from plenum chamber measurements. 

Ideal thrust of total primary exhaust system based on 
measured mass-flow rate, lbs. 

Ideal thrust of total primary exhaust system based on mass- 
flow rate obtained from plenum chamber measurements, lbs. 
Ideal thrust of individual engines (where I = engine number 
(up to 4» based on mass-flow rate obtained from individual 
plenum chamber measurements, lbs. 

Primary exhaust flow air flowmeter frequency, hertz. 
Secondary flow air flowmeter frequency, hertz. 

Ideal thrust based on in-line (not MCV) venturi mass flow, 
where I = engine number, lbs. 

Ratio of specific heats for primary exhaust flow. 

Intercept to be used for determining each engine mass-flow 
rate from plenum chamber measurements, where I = engine 
number. 

Flowmeter number for primary flow air flowmeter. 
Flowmeter number for secondary flow air flowmeter. 
Constant used in chamber mass-flow calculation, used if 
second order curve fit is required, where I = engine number. 


PAGE blank not filmed 
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SYMBOL 

NOMENCLATURE 

KBL 

If set to 1, tertiary flow computation is done. 

If set to 0, tertiary flow computation is omitted. 

KCH(I) 

Slope to be used for determining each engine mass-flow rate 
from plenum chamber measurements, where I = engine 
number. 

KI1 

Internally computed constant. 

KI2 

Internally computed constant. 

KI3 

Internally computed constant. 

KJ1 

Internally computed constant (function of GAMJ). 

KJ2 

Internally computed constant (function of GAMJ). 

KJ3 

Internally computed constant (function of GAMJ). 

KJ4 

Internally computed constant (function of GAMJ). 

KJ5 

Internally computed constant (function of GAMJ). 

KPAV(I) 

Constants used to determine average primary jet total 
pressure ratio from all engines, where I = engine number 
these constants must equal 0.0 or 1.0. 

KPBL(J) 

Constants used to determine average static pressure in 
tertiary duct, where J = probe number. Must equal 0.0 or 


1 . 0 . 


KPCH(I) 

KPS 

KPS(J) 

KPT(I,J) 

KPTBL(J) 


KPTS(J) 


Break pressure for calculation of WPENG(I) for second order 
equations, lbs/sq. in. 

Secondary flowmeter constant (Internally computed). 
Constants used to determine average static pressure in 
secondary air duct, where J = probe number. Must equal 0.0 
or 1.0. 

Constants used in computing jet total pressure, where I = 
engine number and J = probe number. These constants must 
equal 0.0 or 1.0. 

Constants used to determine average total pressure in 

tertiary duct, where J = probe number. Must equal 0.0 or 

1 . 0 . 

Constants used to determine average total pressure in 
secondary air duct, where J = probe number. Must equal 0.0 
or 1.0. 
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SYMBOL 

KR(I,J) 


KSEC 

KSW 


KTAV(I) 

KTT(I,J) 

KV 

KVA(I) 

KVARI(L) 


KVARI(L,I) 


* 


MBLDOT 

MCV 

MDOT 


NOMENCLATURE 

Rake constant for each probe in each engine, where I = 
engine number and J = probe number. If no correction is to 
be made to total pressure probe, then its value should be set 
to 1.0. If probe is bad or does not exist, then its value 
should be set to 0.0. 

If set to 1, secondary flow computation is done. If set to 0, 
secondary flow computation is omitted. 

Switch for chamber, venturi or flowmeter. 

=-l, Venturi mass-flow calculation. 

=0, Flowmeter mass-flow calculation. 

=1, Chamber mass-flow calculation. 

=2, In-line venturis. 

Constants used to determine average primary jet total 
temperature from all engines, where I = engine number. 
These constants must equal 0.0 or 1.0. 

Constants used in determining primary jet total 
temperature, where I = engine number and J = probe 
number. These constants must equal 0.0 or 1.0. 

Venturi constant, used to account for different venturi 
calibrations. It includes venturi throat area and discharge 
coefficient. 

Constants used to determine average static pressure of 
multiple critical venturi, where I = probe number. 

Constants used in the computation of in-line (not MCV) 
venturi weight flow rate, where L = 1 to 4 represents values 
of Pj./P at A/A* of venturi to convert measured static 
pressure at throat to a total pressure and L = 5 to 8 
represents averaging factors (must be 0.0 or 1.0). 

Constants used to associate which in-line (not MCV) venturi 
weight flow rate is related to proper engine, where L = 
venturi number and I = engine number. 

Tertiary mass-flow rate, slugs/sec. 

Venturi meter number. 

Primary mass-flow rate as measured by flowmeter, 
slugs/sec. 
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B-4 

SYMBOL 

MDOTCH 

MSDOT 

NPTE 

NTTE 

NUMENG 

PBL(J) 

PBLAVE 

PCH(I) 

PCHOKE 

PFM 

PFMS 

PS(J) 

PSEC 

PTBL(J) 

PTBLAV 

PTB/PTJ 

PTB/PTO 

PTENG(I) 

PTENG(I)/PO 

PTENGO(I) 


NOMENCLATURE 

Primary mass-flow rate as computed from plenum chamber 
measurements, slugs/sec. 

Secondary flow mass-flow rate, slugs/sec. 

Number of total pressure probes in each engine, where I = 
engine number. (Internally computed). 

Number of total temperature probes in each engine, where I 
= engine number. (Internally computed). 

Number of engines in model (maximum of 4). NUMENG = 0 
for aerodynamics model (no other constants required). 

Static pressure measurements in the tertiary duct (up to 4), 
where J = probe number, lbs/sq. in. 

Average static pressure in the tertiary duct, lbs/sq. in. 
Individual engine-plenum-chamber total pressure, I = engine 
number, lbs/sq. in. 

Primary jet-total-pressure ratio for choked flow. 

Pressure measured at primary flow flowmeter, lbs/sq. in. 
Pressure measured at secondary flow flowmeter, lbs/sq. in. 
Static pressure measurements in the secondary flow duct (up 
to 4), where J = probe number, lbs/sq. in. 

Average static pressure in the secondary flow duct, lbs/sq. 
in. 

Total pressure measurements in the tertiary duct (up to 4), 
where J = probe number, lbs/sq. in. 

Average total pressure in the tertiary duct, lbs/sq. in. 

Ratio of tertiary total pressure to primary jet total 
pressure. 

Ratio of tertiary total pressure to free-stream total 
pressure. 

Average primary jet total pressure in each engine, where I = 
engine number, lbs/sq. in. 

Ratio of average primary jet total pressure in each engine 
to tunnel static pressure, where I = engine number. 

Ratio of average primary jet total pressure in each engine 
to tunnel static pressure, where I = engine number. 
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SYMBOL 

PTJ(I,J) 

PTJ/PO 

PTS(J) 

PTS/PTJ 

PTS/PTO 

PTSEC 

PTV 

PV 

PV1 

PV2 

PV/PTV 

PVEN(I) 


PVRI(L) 

RJ 

RS 

RV 

TCH(I) 

TFM 

TFMS 

THETBL 

THETSE 

TTBL 


NOMENCLATURE 

Individual primary jet total pressure measurements, where I 
= engine number and J = probe number, lbs/sq. in. 

Average primary jet total pressure ratio (all engines). 

Total pressure measurements in the secondary flow duct, 
where J = probe number, lbs/sq. in. 

Ratio of secondary flow total pressure to primary jet total 
pressure. 

Ratio of secondary flow total pressure to free-stream total 
pressure. 

Average total pressure in the secondary flow duct, lbs/sq. 
in. 

Tertiary venturi total pressure, lbs/sq. in. 

Tertiary venturi static pressure, lbs/sq. in. 

Averaged multiple critical venturi static pressure upstream 
of venturi throat, lbs/sq. in. 

Averaged multiple critical venturi static pressure 
downstream of venturi throat, lbs/sq. in. 

Ratio of tertiary venturi static pressure to tertiary total 
pressure. 

Multiple critical static pressure, where 1=1 and 3 are 
upstream and 1 = 2 and 4 are downstream of venturi throat, 
lbs/sq. in. 

In-line (not MCV) venturi static pressure, where L = venturi 
number, lbs/sq. in. 

Gas constant for primary flow, ft/degree Rankine. 

Gas constant for secondary flow, ft/degree Rankine. 

Gas constant for tertiary flow, ft/degree Rankine. 

Individual engine-plenum chamber total temperature, I = 
engine number, °F. 

Temperature at primary flowmeter, °F. 

Temperature at secondary flowmeter, °F. 

Tertiary flow corrected mass-flow ratio. 

Secondary flow corrected mass-flow ratio. 

Total temperature of tertiary flow, °F. 


19 



B-6 

SYMBOL 

TTENG(I) 

TTJ(I,J) 

TTJAVG 

TTSEC 

TTV 

TV 

TVRI(L) 

VRATIO 

WI 

WIENG(I) 

WMCV 

WMCV/WI 

WP 

WPBL 

WPCHR 

WPCHR/WI 


WPENG(I) 


WPSEC 

WP/WI 

WPE/WIE(I) 


NOMENCLATURE 

Average primary jet total temperature in each engine where 
I = engine number, °F. 

Individual primary jet total temperature measurements 
where I = engine number and J = probe number, °F. 

Average primary jet total temperature (all engines), °F. 
Secondary flow total temperature, °F. 

Temperature at the tertiary venturi, °F. 

Multiple critical venturi temperature, °F. 

Temperature at the in-line (not MCV) venturi, where L = 
venturi number, °F. 

Ratio of multiple critical venturi static pressures (should be 
less than 0.93). 

Ideal weight flow of primary flow, lbs/sec. 

Ideal weight flow of each individual engine primary flow, 
where I = engine number, lbs/sec. 

Multiple critical venturi weight flow rate, lbs/sec. 

Ratio of multiple critical venturi weight flow rate to ideal 
weight flow rate. 

Measured weight flow of air primary flow flowmeter or 
venturi, lbs/sec. 

Tertiary weight flow rate obtained from venturi, lbs/sec. 
Total primary flow weight flow rate obtained from plenum 
chamber measurements, lbs/sec. 

Discharge coefficient of total primary flow system as 
obtained from plenum chamber measurements for entire 
system. 

Primary flow weight flow rate of each engine obtained from 
plenum-chamber measurements, where I = engine number, 
lbs/sec. 

Secondary flow weight flow rate, lbs/sec. 

Primary flow discharge coefficient using flowmeter or 
venturi weight flow rate for entire system. 

Discharge coefficient of each individual engine as obtained 
from plenum-chamber measurements, where I = engine 
number. 
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SYMBOL 

WPVRI 

WPVRI/WI 

WV/WKI) 

WVRI(L) 

Z 

ZS 


NOMENCLATURE 

Sum of in-line (not MCV) venturi weight flow rates, lbs/sec. 
Ratio of summation of in-line (not MCV) venturi weight flow 
rate to ideal weight flow rate. 

Ratio of in-line (not MCV) venturi weight flow to ideal 
weight flow of each engine, where I = engine number. 

In-line (not MCV) venturi weight flow rate, where L = 
venturi number, lbs/sec. 

Primary flowmeter constant. (Internally computed). 
Secondary flowmeter constant. (Internally computed). 
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APPENDIX B 
Module B 

Jet Exhaust Measurements 

A. Required Constants 

1. All constants are initialized to a value of zero. The project engineer 
needs to supply only those constants which are required for the quantities 
to be computed. In addition, by logical use of combinations of these 
constants, several options are available to the project engineer. One of 
these options is discussed later. 

2. NUMENG - number of engines in model. NUMENG = 0 for aerodynamics 
model (no other constants are required). 

3. KR(I,J) - Rake constant for each probe in each engine, where I = engine 

number and J = probe number. 

If no correction is to be made to the total pressure probe, then its value is 
set equal to 1.0. If the probe is faulty or does not exist, then its value is 
set equal to 0.0. 

Example: Two engines; five probes in the first, and three probes in the 

second. 

Engine 1 is corrected to integrated rake values, engine 2 
probes are uncorrected. 

KR(1,1) = 1.051 
KR(1,2) = .986 
KR(1,3) = .972 
KR(1,4) = .987 
KR(1,5) = 1.058 
KR(2,1) = 1.0 
KR(2,2) = 1.0 
KR(2,3) = 1.0 

nwcrorc PAGE sum WT «J*° 
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Note that there is no need to supply those constants which equal zero since they are 
assumed to be zero if not supplied. 

4. Special Case: A twin-engine configuration with only one set of 

chamber measurements is not uncommon. The following constants 
are used. 

NUMENG = 2 

AENG(l) = total orifice nozzle area 
AENG(2) = 0.0 

This combination of constants yields the following, nonstandard, results: 

WPENG(l) = total weight flow based on pressure and temperature 
measurements of engine 1. 

WPENG(2) = 0.0 

WPE/WIE(1) and WPE/WIE(2) are meaningless 

FIENG(l) = total ideal thrust based on pressure and temperature 
measurements of engine 1. 

FIENG(2) = 0.0 

WPCHR, MDOTCH, WPCHR/WI, FICHR, and CFICHR are based on 
pressure and temperature measurements in engine 1 rather than on 
the average values of both engines. 

B. Test for Exhaust Model 

1. The constant required from the project engineer is NUMENG (0 to 4). 

IF NUMENG = 0, skip module B. 

C. Compute Common Constants 

1. The constants required from the project engineer are GAMJ and RJ. 



B— 1 1 


KJ1 = 


( t—\ 

\GAMJ + 1/ 


GAMJ + 1 
2(GAMJ - 1) 


a /GAMJ^J2a74 
V RJ 


(Eq. B-l) 


GAMJ * 64.348 
(GAMJ - 1)RJ 


(Eq. B-2) 


, / 2 * (GAMJ) * (RJ) 

V (GAMJ - 1) * 32.174 


(Eq. B-3) 


KJ4 = 


GAMJ - 1 
GAMJ 


(Eq. B-4) 


KJ5 = 


1 _ 

GAMJ 


(Eq. B-5) 


GAMJ 

r /gamj - ~ i 

PCHOKE = 1 + ^ M 


(Eq. B-6) 


D. Individual Engine Measurements 

1. This permits computation for four separate engines with the following 
instrumentation in each engine: 

a. jet total pressures 

b. jet total temperatures 

c. chamber pressure 

d. chamber temperature 

2. Jet total pressure 

a. Jet total pressure will always be called PTJ(I,J), where I = 
engine number and J = probe number. 
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b. An example of representing the third measurement (probe 3) 
of jet total pressure in engine 2 is named PTJ(2.3). 

c. The constants required from the project engineer are KR(I,J) 
and KPT(I,J). 

NPTE(I) 

l PTJ(I,J) * KR(I,J) 

pteno(i) = fffygny (Eq - B ' 7> 

l KPT(I,J) 

J = 1 

PTENGO(I) = PT pQ G --- (Eq. B-8) 

3. Jet total temperature 


a. Jet total temperature measurements are always called 
TTJ(I,J), where I = engine number and J = probe number. 


b. An example of the first measurement (probe 1) of jet total 
temperature in engine 3 is named TTJ(3.1). 


c. The constants required from the project engineer are KTT(I,J) 
and NTTE(I). 


NTTE(I) 

l TTJ(I,J) * KTT(I,J) 

TTENG(I) = UjSgfiT 

l KTT(I,J) 

J = 1 


(Eq. B-9) 


4. Chamber weight flow for each engine. 

a. The constants required from the project engineer are KCH(I), 
ICH(I), KAE(I), AT(I), AENG(I) and KPCH(I). 
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If PCH(I) < KPCH(I) 
then 

WPENG(I) = AEN G(I) * PCH(I) * KJ1 *[ ICH(I) + KCH(I) * PCH(I) + KAE(I) * PCH(I) 2 

yTCH(I) + 459.67 

If PCH(I)> KPCH(I) then 


WPENG(I) = 


_ AENG(I) * PCH(I) * KJ1 *[ ICH(I + 4) + KCH(I + 4) * PCH(I) + KAE(I + 4) * PCH(I)0 


n/tCH(I) 


+ 459.67 


(Eq. B-10) 


5. Ideal weight flow for each engine. 


a. The nozzle choke total pressure ratio is calculated internally 
and is called PCHOKE. 

b. The constant required from the project engineer is AT(I). 

If PTENGO(I) is greater than PCHOKE, use equation B-ll. 


WIENG(I) = 


CkjJ * CpTENGdll * CaT(I)) 
yTTENG(I) + 459.67 


(Eq. B-ll) 


If PTENGO(I) is less than or equal to PCHOKE, use equation B-12. 


KI1 = 


KJ2 

(TTENG(I) + 459.67) 



1 \ KJ4 

PTENGO(I) / 


If KI1 is less than 0, KI1 = .0001 


(Eq. B-12) 
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then 


WIENG(I) = n/kU * AT(I) * PTENG(I) * (pflifew) KJ5 (Eq ‘ B " 13) 


Note to the project engineer: If the engine is shrouded, then a local static pressure in 
the nozzle shroud should be used rather than PO. The engineer must supply a new 
equation for KI1 and WIENG(I). 


6. Discharge coefficient for each engine based on chamber weight flow. 


WPE/WIE(I) = 


WPENG(I) 

WIENG(I) 


(Eq. B-14) 


If WIENG(I) = 0, WPE/WIE(I) = 0 


7. Ideal thrust for each engine based on chamber weight flow. 


K12 = TTENG(I) + 459 


.67J * [j " 


If K12 is less than 0, KI2 = .0001 


1 KJ4 “| 

PTENGO(I) | 

(Eq. B-15) 


FIENG(I) = . KJ3] * jwPENG(I)J * [7 k12 ] 


(Eq. B-16) 


8. In-line venturi: weight flow for each engine. The equations given below 
are for critical flow venturi and are intended to be very general. 

A(I) = j [vKRI(I,4) * (TVRI(L) + 459.67) + VKRI(I,3)J * 

(TVRI(L) + 459.67) + VKRI(I,2)J * (TVRI(L) + 459.67) + VKRI(I,1) 

A(I) where I = 1 to 4 are constants which go into the compressibility term, 
C*. As seen, a 3rd order equation capability exists. Values of VKRI(I,1) 
to VKRI(I,4) can be input using 'T' cards to allow use of most any critical 
venturi. 
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C* = [(A(4) * PVRI(L) + A(3)) * PVRI(L) + A(2)] * PVRI(L) + A(l) 

TS = 0.8333 * TVRI(L) + 459.67 

VIS = 6.086248 * 10 -8 * (TS) 1,5 / [TS + 198.6) 


Individual venturi mass flow is then computed using 


WVRI(L) = PVRI(L) * KVARI(L) * AVRKL) * g * C* * CDSI(L) 
/g * RJ * (TVRI(L) + 459.67) 


NOTE : CDSI(L) represents the discharge coefficient of individual 

venturi. It is obtained using an iterative scheme based on 
venturi throat Reynolds number. A table of CD versus 
RDUCT is required for each venturi. RDUCT is computed 
using 


RDUCT(L) = WVRI(L)/(AVRI(L) * VIS) 


Because of the complexity of this computation, an example is included. 
The following information is contained within the data reduction program 
when using the twin critical venturis which measure total mass flow in the 
groundstand (B1234). 


VKRI(1,4) = 0.0 
VKRI(1,3) = -1.43545E-8 
VKRI(1,2) = 1.36243E-5 
VKRI(1,1) = 0.68166 
VKRI(2,4) = 0.0 
VKRI(2,3) = 4.49456E-10 
VKRI(2,2) = -6.06496E-7 
VKRI(2,1) = 2.14835E-4 


VKRI(3,4) = 0.0 
VKRI(3,3) = 1.64438E-13 
VKRI(3,2) = -1.90568E-10 
VKRI(3,1) = 5.4424E-8 
VKRI(4,4) = 0.0 
VKRI(4,3) = 0.0 
VKRI(4,2) = 0.0 
VKRI(4,1) = 0.0 
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KVARI(l) = 1.0040 AVRI(l) = .272009 KVARI(5) = 1.0 

KVARI(2) = 1.0039 AVRI(2) = .264481 KVARI(6) = 1.0 


Only the KVARI and AVRI constants are required to be input by an 
engineer. Both venturis use the same CD versus RDUCT relationship, 
which is not a table lookup but simply a second order equation. Of course 
a table lookup could be used in lieu of the equation. 

The CDSI equation for twin critical venturis in groundstand: 

CDSI(L) = 0.993507 + 3.5062E-4(RDUCT(L)) - 1.1269E-5(RDUCT(L) 2 ) 

where RDUCT(L) = WVRI(L)/(AVRI(L) * VIS * 1.0E06) 

9. Discharge coefficient for each engine based on in-line venturi weight 
flow. 

WV/WIU) = W VRI(I)/ WIEN G(I) 

10. Ideal thrust for each engine based on in-line venturi weight flow. 

FVRI(I) = WVRI(I) * KJ3 * J K12 (Eq. B-17) 

Total Exhaust System Properties 
1. Average total pressure ratio. 

a. The constant required from the project engineer is KPAV(I). 


NUMENG 

l [KPAV(I) * PTENGO(I)] 

PWPO. (Eq. B-18) 

l KPAV(I) 

1 = 1 



2. Average total temperature. 
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a. The constant required from the project engineer is KTAV(I). 


TTJAVG = 


NUMENG 

l (kTAV(I) * TTENG(I)} 

I = 1 

NUMENG 
I KTAV(I) 

1 = 1 


3. Total weight or mass flow. 


(Eq. B-19) 


a. The total system weight flow is in units of lb/sec. 

b. The total system mass flow is in units of slugs/sec. 

c. The constants required from the project engineer are: 

(1) INTFM1 and MCV 

(2) KSW selects mass flow computation 

= 1; chamber flow 
= 0; flowmeter 
=-l; MCV venturi 
= 2; in-line venturis 

If KSW = 1 (chamber mass flow calculation) 

(Eq. B-20) 
(Eq. B-21) 


NUMENG 

WPCHR = l WPENG(I) 
1 = 1 


MDOTCH = 


WPCHR 

32.174 
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If KSW = 0 (air model with flowmeter) 


Z and KP are determined from standardized flowmeter tables 


(FM1) * (PPM) * (144. ) 

WP = (RJ) * (Z) * (KP) * (TFM + 459.67) 


(Eq. B-22) 


MDOT = 


WP 

32.174 


(Eq. B-23) 


If KSW = -1 (venturi mass flow calculation), the venturi code, MCV, is decoded 
to derive those venturi present 


PV1 = 


KVA1 * PVEN1 + KVA3 * PVEN3 
KVA1 + KVA3 ” 


PV2 = 


KVA2 * PVEN2 + KVA4 * PVEN4 
KVA2 + KVA4 


PV2 

VRATIO = pyj 


A(I) = ((VK(I,4) * TV + VK(I,3)) * TV + VK(I,2)) * TV + VK(I,1) (Eq. B-24) 

C* = ((A(4) * PV1 + A(3)) * PV1 + A(2)) * PV1 + A(l) (Eq. B-25) 

TS = 0.8333 * (TV + 459.67) ( E< 3* B-26 ^ 

VIS = 6.086248 * 10' 8 * (TS) 1,5 /(TS + 198.6) (Eq. B-27) 

\ 1 /2 

WMCV = l PV1 * AREAV(I) * (C*) *^V+ 4 459.67)Rj) * CD(I) 

I (Eq. B-28) 

ARMCV = l AREAV(I) ( E ^ B “ 29 ^ 

I 

The above summations are over the venturi present. CD(I) is 
computed by linear interpolation from a table of CD vs 
RNMCV 


where 

RNMCV = WMCV/(ARMCV*VIS) 


(Eq. B-30) 


B-19 

An iterative scheme is used until successive computations of 
WMCV differ by a desired accuracy. 

4. If KSW = 2 (in-line venturis) 

4 

WPVRI = £ WVRI(L) * KVARI(L + 4) (Eq. B-31) 

L = 1 

5. Ideal weight flow (total). 

a. Ideal weight flow of the total system is computed 
NUMENG 

WI= l WIENG(I) (Eq. B-32) 

1 = 1 

6. Discharge coefficient for the entire system. 

a. The discharge coefficient using weight flow from a flowmeter 
or a venturi and the discharge coefficient using weight flow 
from chamber measurements are computed. 


(Eq. B-33) 


(Eq. B-34) 
(Eq. B-35) 


If KSW = 2 WP = WPVRI 
KSW = 1 WP = WPCHR 
KSW = 0 WP = WP 
KSW = -1 WP = WMCV 


MDOT = 


WP 


32.174 


WP/WI = 


WP 

WI 


WPCHR/WI = 


WPCHR 

WI 
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WMCV/WI = 
WPVRI/WI = 


WMCV 

WI 

WPVRI 

WI 


If WI = 0; WP/WI = WPCHR/WI = WMCV/WI = WPVRI/WI = 0 

(Eq. B-36) 


Ideal thrust for the entire system. 


a. The ideal thrust, FICHR, and ideal thrust coefficient CFICHR 
are obtained from chamber weight flow. 


b. The ideal thrust, FI, and ideal thrust coefficient CFI are 
obtained from flowmeter or venturi measured weight flow. 

c. Note that MACH, PO and QO are from Module A. 


d. The constant required from project engineer is AREF 


NUMENG 

FICHR = l FIENG(I) 

I = 1 

If MACH< .1, 


CFICHR = 


FICHR 

(PO) * (AREF) 


(Eq. B-37) 


(Eq. B-38) 


If MACH> .1, 


CFICHR = 


FICHR 


(QO) * (AREF) B - 39 > 


KI3 = (TTJAVG + 459.67) * 


1- — L 


KJ4 


PTJ/PO 


(Eq. B-40) 


If KI3 < 0; KI3 = .0001 
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FI = (KJ3) * (WP) * (VKf3 ) (Eq. B-41) 

WP from flowmeter if KSW = 0 
WP from venturi if KSW = -1 
WP = WPCHR if KSW = 1 
WP = WPVRI if KSW = 2 

If MACH < .1, 


CFI (PO) * (AREF) (Eq - B_42) 


If MACH > .1, 

CF1= ^ V (ARE F) (Eq. B— 43) 

If KSEC = 0, skip equations B-44 through B-50. 

F. Secondary Flow Measurements 

1. Secondary passage total pressure. 

a. The total pressure measurements PTS(J) in the secondary air 
passage (up to 4) are used to compute the average secondary 
passage total pressure. 

b. The constant required from the project engineer is KPTS(J). 

4 

l KPTS(J) * PTS(J) 

PTSEC = — (Eq. B-44) 

l KPTS(J) 

J = 1 
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2. Secondary passage static pressure. 

a. Static pressure measurements PS(J) in the secondary air 
passage (up to 4) are used to compute the average static 
pressure in the secondary air passage. 

b. The constant required from the project engineer is KPS(J) 

4 

l KPS(J) * PS(J) 

PSEC = ^ (Eq. B-45) 

l KPS(J) 

J = 1 


3. Secondary duct total temperature. 


a. The total temperature TTSEC in the secondary duct is handled 
in the standard program for quantities. 


4. Secondary mass flow. 


a. The constants required from the project engineer are RS, KPS, 
ZS, INTFMS. KPS and ZS are determined internally from 
INTFMS constant. 


WPSEC = 


(FMS) * (PFMS) * (144.0) 

(RS) * (ZS) * (KPS) * (TFMS + 459.67) 


, lbs/sec 


(Eq. B-46) 


MSDOT = 


f£f£r ’ slugs/sec 


(Eq. B-47) 


5. 


Pumping characteristics 


PTS/PTJ = 


PTSEC 

(PTJ/PO) * (PO) 


(Eq. B-48) 
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PTS/PTO = 


PTSEG 

PTO 


(Eq. B-49) 


6. Corrected mass flow ratio 


THETSE = 


MSDOT / (TTSEC + 459.67) * RS 

MDOT V (TTJAVG + 459.67) * RJ 


(Eq. B-50) 


If KBL = 0, skip equations B-51 through B-57. 


G. Tertiary Flow Measurements 

1. Tertiary duct total pressure. 


a. The total pressure measurements PBL(J) in the tertiary duct 
(up to 4) are used to compute the average tertiary duct total 
pressure. 

b. The constant required from the project engineer is KPTBL(J). 

4 

J KPTBL(J) * PTBL(J) 

PTBLAV = (Eq. B-51) 

J KPTBL(J) 

J = 1 


2. Tertiary duct static pressure. 

a. Static pressure measurements PBL(J) in the tertiary duct (up 
to 4) are used to compute the average static pressure in the 
tertiary duct. 
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b. 


The constant required from the project engineer is KPBL(J). 


4 

y KPBL(J) * PBL(J) 

PBLAVE = — ( E <1* B " 52 ) 

y KPBL(J) 

J = 1 


3. Tertiary duct total temperature. 

a. Total temperature in the tertiary duct TTBL is handled in the 
standard program for quantities. 


4. Tertiary mass flow. 

a. Venturi total pressure, PTV, and venturi static pressure, PV, 
are required. 

b. Tertiary weight flow is in units of lbs/sec. 


Tertiary mass flow is in units of slugs/sec. 


d. The constants required from the project engineer are RV, KV. 


PV/PTV = 


PV 

PTV 


WPBL = KV 


/ P V W PTV \ 
\Po a o ] \yTTV + 459.67 / 


(Eq. B-53) 


where ■ ft — is a function of PV/PTV and is 
Po a o 

determined from slopes and intercepts supplied by 
transonic tunnel personnel. 


the 16-foot 
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MBLDOT - , slugs/sec (Eq. B-54) 


5. Pumping characteristics. 


PTB/PTJ = 


PTBLAV 

(PTJ/PO) 


(Eq. B-55) 


PTB/PTO = 


PTBLAV 

PTO 


(Eq. B-56) 


6. Corrected mass flow ratio. 


THETBL = 


MBLDOT 

MDOT 



(TTBL + 459.67) * RV 
(TTJAVG + 459.67) * RJ 


(Eq. B-57) 
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MODULE C SKIN FRICTION DRAG 
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SYMBOL 

AREF 

AWET(I) 


CDF 

CDF AFT 

CDFF 

CDFHT 

CDFNOZ 

CDFR(I) 

CDFVT 

CDFW 

FL(I) 


FORMF(I) 


NOMENCLATURE 

Model reference area used for coefficients, sq. in. If 
module B is used, this constant is already specified. 

Model wetted areas, sq. in. 

Where AWET(l) = total fuselage wetted area. 

AWET(2) = fuselage wetted area up to metric break. 
AWET(3) = fuselage wetted area up to nozzle connect 
station. 

AWET(4) = wing wetted area. 

AWET(5) = vertical tail wetted area. 

AWET(6) = horizontal tail wetted area. 

AWET(7) = optional, for additional body. 

Total skin friction drag coefficient. 

Afterbody plus nozzle skin friction drag coefficient. 

Total fuselage skin friction drag coefficient. 

Horizontal tails (canards) skin friction drag coefficient. 
Nozzle skin friction drag coefficient. 

Individual skin friction drag coefficients calculations. 
Vertical tails(s) skin friction drag coefficient. 

Wing skin friction drag coefficient. 

Model reference lengths, feet. 

Where FL(1) = fuselage length. 

FL(2) = fuselage length up to metric break. 

FL(3) = fuselage length up to nozzle connect station. 
FL(4) = wing mean aerodynamic chord. 

FL(5) = vertical tail mean aerodynamic chord. 

FL(6) = horizontal tail mean aerodynamic chord. 

FL(7) = optional 
Form factors 

Where FORMF(l) = fuselage. 

FORMF(2) = wing. 

FORMF{3) = vertical tail 
FORMF(4) = horizontal tail. 

FORMF(5) = optional. 
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SYMBOL 

NOMENCLATURE 

KFAFT 

Constant used to include proper terms in total skin friction 


drag term, CDF. Must equal 0.0 or 1.0. 

KFF 

See KFAFT. 

KFNOZ 

See KFAFT. 


r- u 
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APPENDIX C 
Module C 

Skin Friction Drag 

Skin friction drag is computed by the method of Frankl and Voishel 1 for compressible, 
turbulent flow on a flat plate. 

A. Required Constants 

All constants are initialized to a value of 0.0 except FORMF(I) which is 
initialized to a value of 1.0. 

1. AWET(I) 

2. FORMF(I) 

Form factors may be obtained from LWP - 1120. 

Fuselages FORMF(l) = 1.0 + 1.5(d/l) 1,5 + 7(d/l) 3 (Eq. C-l) 

Empennages FORMF(I) = 1.0 + 1.44(t/c) + 2(t/c) 2 (Eq. 02) 

3. The model reference lengths (FL(I)), are given in the nomenclature 
section. 

.2 

4. The model reference area (AREF) is used for coefficients, in . If jet 
exhaust measurements are used, this constant is already specified. 

5. The constants (KFF, KFAFT, KFNOZ) used to include proper terms in 
total skin friction drag term, CDF, must equal 0 or 1. 

B. Test for Skin Friction Calculation 

If AWET(l) = 0, skip the calculations for the skin friction drag in this module. 


1 Frank!, F., and Voishel, V. Friction in the turbulent boundary layer of a compressible gas 

at high speeds. TM NACA No. 1032, 1942. 
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C. Fuselage Skin Friction Drag 

1. The constants required from the project engineer are AWET(l), AWET(2), 
AWET(3), FL(1), FL(2), FL(3), AREF, and FORMF(l). 

J= 3 


If AWET(2) = 0 and AWET(3) = 0, J = 1 
If AWET(2) 4 0 and AWET(3) = 0, J = 2 


Calculate CDFR(I) for I = 1, J 


CDFR(I) = 2 4fi7 7 

(1 + .2 MACHT *j lo g 10 

(RN/FT) * FL(I) 

) 2.58 

> * AREF 

\ 


' (Eq. C-3) 

If MACH < .1, CDFR(I) = 0.0 



CDFF = CDFR(l) 


(Eq. C-4) 

If AWET(2) 4 0, 



CDF AFT = CDFR(l) - CDFR(2) 


(Eq. C-5) 

If AWET(3) 4 0, 



CDFNOZ = CDFR(l) - CDFR(3) 


(Eq. C-6) 


D. Empennage Skin Friction Drag 

1. The constants required from the project engineer are AWET(4), AWET(5), 
AWET(6), FL(4), FL(5), FL(6), AREF, FORMF(2), FORMF(3), FORMF(4), 
KFF, KFAFT, and KFNOZ. 
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Calculate CDFR(I) for I = 4, 7 


J = 1-2 


CDFR(I) = 


.472 » AWET(I) « FORMF(J) 


(1 


+ .2 MACH 2 )* 467 


*| lo «10 


(RN/FT) * FL(I) 


2.58 

* AREF 
(Eq. C-7) 


IF MACH < . 1 , CDFR(I) = 0 
CDFW = CDFR(4) 


CDFVT = CDFR(5) 
CDFHT = CDFR(6) 


Total Skin Friction Drag 

1. Skin friction drag of the entire model is computed. 

CDF = (KFF * CDFF) + (KFAFT * CDFAFT) + (KFNOZ * CDFNOZ) 

+ CDFW + CDFVT + CDFHT + CDFR(7) (Eq. C-8) 
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MODULE D BALANCE LOADS AND MODEL ATTITUDES 


D-l 


SYMBOL NOMENCLATURE 

The arrays FO through F20 are forces and moments whose units are lbs and in. Ibs. 


AF(I,J) 

AFO(I) 

AFT(I) 

AFTARE(I) 

ALPHA 

AMOM(I) 

ARB(H,K) 


ARP(n,K) 


ARPB(II,K) 



Axial force, lbs., where I = balance number and J = 
correction number. 

Initial axial load, lbs., where I = balance number. 

Total axial load, lbs., where I = balance number. 

Axial weight tares, lbs., where I = balance number. 

Model angle of attack, degrees. 

Axial force momentum correction, lbs., where I = balance 
number. 

Areas or momentum arms * areas used with PBASE(H) for 
computing base force and base moment tares. Care should 
be used to insure proper tare force signs. Area and arm 

units must be consistent with units of base pressures and 

balance components. Second balance, sq. in., where K = 
component number and II = orifice number. 

Areas or momentum arms * areas used with PBASE(II) for 
computing base force and base moment tares. Care should 
be used to insure proper tare force signs. Area and arm 

units must be consistent with units of base pressures and 

balance components. Third balance, sq. in., where K = 
component number and II = orifice number. 

Areas or momentum arms * areas used with PBASE(II) for 
computing base force and base moment tares. Care should 
be used to insure proper tare force signs. Area and arm 
units must be consistent with units of base pressures and 
balance components. First balance, sq. in., where K = 
component number and II = orifice number. 

Initial balance loads, axial force, lbs. (Weight Tares) 

Balance component quantity corrected for high interactions 
coupled with high model restraints, axial force, lbs. (Weight 
Tares) 


PRECEDING PAG€ BLAP<k NOT FILMED 
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SYMBOL 

A 4 

BETA 

BSPAN(I) 

CA(I) 

CABASE(I) 

CAREF(I) 

CC(I) 

CD(I) 

CDBASE(I) 

CDS(I) 

CHORD(I) 

CL(I) 

CLS(I) 

CLSQR(I) 

CMX(I) 

CMXREF(I) 

CMXS(I) 

CMXW(I) 

CMY(I) 

CMYREF(I) 


NOMENCLATURE 

Balance component quantities corrected for balance 
orientation to gravity axis, axial force, lbs. (Weight Tares) 
Angle of sideslip, degrees. 

Roll and yaw moments reference length. Normally wing 
span, inches, where I = balance number. 

Axial force coefficient in the body axis, where I = balance 
number. 

Base axial force coefficient, where I = balance number. 

Axial force coefficient in the reference axis, where I = 
balance number. 

Crosswind coefficient in the wind axis, where I = balance 
number. 

Drag coefficient in the wind axis, where I = balance number. 
Base drag coefficient, where I = balance number. 

Drag coefficient in the stability axis, where I = balance 
number. 

Pitching moment reference length. Normally wing mean 
aerodynamic chord, inches, where I = balance number. 

Lift coefficient in the wind axis, where I = balance number. 
Lift coefficient in the stability axis, where I = balance 
number. 

Lift coefficient squared, where I = balance number. 

Rolling moment coefficient in the body axis, where I = 
balance number. 

Rolling moment coefficient in the reference axis, where I = 
balance number. 

Rolling moment coefficient in the stability axis, where I = 
balance number. 

Rolling moment coefficient in the wind axis, where I = 
balance number. 

Pitching moment coefficient in the body axis, where I = 
balance number. 

Pitching moment coefficient in the reference axis, where I = 
balance number. 
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SYMBOL 

CMYS(I) 

CMYW(I) 

CMZ(I) 

CMZREF(I) 

CMZS(I) 

CMZW(I) 

CN(I) 

CNBASE(I) 

CNREF(I) 

CPBASE(H) 

CPMBASE(I) 

CRMBASE(I) 

CYBASE(I) 

CY(I) 

CYMBASE(I) 

CYREF(I) 

CYSd) 

Cl 

C2 

AA 


NOMENCLATURE 

Pitching moment coefficient in the stability axis, where I = 
balance number. 

Pitching moment coefficient in the wind axis, where I = 
balance number. 

Yawing moment coefficient in the body axis, where I = 
balance number. 

Rolling moment coefficient in the reference axis, where I = 
balance number. 

Yawing moment coefficient in the stability axis, where I = 
balance number. 

Yawing moment coefficient in the wind axis, where I = 
balance number. 

Normal force coefficient in the body axis, where I = balance 
number. 

Base normal force coefficient, where I = balance number. 
Normal force coefficient in the reference axis, where I = 
balance number. 

Base pressure coefficient, where n = orifice number. 

Base pitching moment coefficient, where I = balance 
number. 

Base rolling moment coefficient, where I = balance number. 
Base side force coefficient, where I = balance number. 

Side force coefficient in the body axis, where I = balance 
number. 

Base yawing moment coefficient, where I = balance number. 
Side force coefficient in the reference axis, where I = 
balance number. 

Side force coefficient in the wind axis, where I = balance 
number. 

Linear balance interactions. 

Nonlinear balance interactions. 

W(AF), axial force weight tares, lbs. 

WY(RM), rolling moment weight tares, in. lb. 
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SYMBOL 

“l 
A rrij 

Am. 

an 


DPBASE(n) 
A W(I) 

Ay 

FA 

FA(I) 

FA(I,L) 

FABASE(I) 

FAMAX 

FAMOM(I) 

FAREF'(I) 

FAREF(I) 

FC(I) 

FD(I) 

FDS(I) 

FL(I) 

FLS(I) 

FN 

FNBASE(I) 

FN(I) 


NOMENCLATURE 

WZ(RM), polling moment weight tares, in. lb. 

WX(PM), pitching moment weight tares, in. lb. 

WZ(PM), pitching moment weight tares, in. lb. 

W(NF), normal force weight tares, lbs. 

WX(YM), yawing moment weight tares, in. lb. 

WY(YM), yawing moment weight tares, in. lb. 

Differential base pressures, where II = orifice number. 

Half weight of balance, lbs., where I = balance number. 
Used in weight tares program. 

W(SF), side force weight tares, lbs. 

Axial force, lb. 

Final body axis axial force, lbs., where I = balance number. 
Balance axial force rotated (L = 1) and translated (L = 2) to 
body axis, where I = balance number. 

Base axial force, lbs., where I = balance number. 

Maximum absolute value of axial force, lbs. 

Axial force due to momentum of flow, lbs., where I = 
balance number. 

Axial force rotated to reference axis, lbs., where I = balance 
number. 

Axial force translated to reference axis, lbs., where I = 
balance number. 

Crosswind force in the wind axis, lbs., where I = balance 
number. 

Drag force in the wind axis, lbs., where I = balance number. 
Drag force in the stability axis, lbs., where I = balance 
number. 

Lift force in the wind axis, lbs., where I = balance number. 
Lift force in the stability axis, lbs., where I = balance 
number. 

Normal force, lb. 

Base normal force, lbs., where I = balance number. 

Final body axis normal force, lbs., where I = balance 
number. 
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SYMBOL 

FN(I,L) 

FNMAX 

FNREF'(I) 

FNREF(I) 

FP 

FT 

FTARE 

FUT 

FY 

FY(I) 

FY(I,L) 

FYBASE(I) 

FYMAX 

FYREF'(I) 

FYREF(I) 

FYS(I) 

FO 

FI 

F2 

F3 

F4 


F5 


F6 


NOMENCLATURE 

Balance normal force rotated (L = 1) and translated (L = 2) 
to body axis, where I = balance number. 

Maximum absolute value of normal force, lbs. 

Normal force rotated to reference axis, lbs., where I = 
balance number. 

Normal force translated to reference axis, lbs., where I = 
balance number. 

All product combinations of vector FT. 

Corrected total loads. 

Tare loads. 

Uncorrected total loads. 

Side force, lbs. 

Final body axis side force, lbs., where I = balance number. 
Balance side force rotated (L = 1) and translated (L = 2) to 
body axis, where I = balance number. 

Base side force, lbs., where I = balance number. 

Maximum absolute value of side force, lbs. 

Side force rotated to reference axis, lbs., where I = balance 
number. 

Side force translated to reference axis, lbs., where I = 
balance number. 

Side force in the stability axis, lbs., where I = balance 
number. 

Initial loads. 

Uncorrected balance quantities. 

Balance component quantities corrected for interactions. 
Vector representing balance component quantities corrected 
for high interactions coupled with high model restraints. 
Vector representing balance quantities corrected for 
balance orientation to gravity axis, attitude loads, and 
weight tares. 

Vector representing balance quantities corrected for method 
of attachment. 

Balance components rotated to the model (body) axis. 
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SYMBOL 

F7 

F8 

F9 

F10 

Fll 

F12 

F13 

F14 

F15 

F16 

F17 

F18 

F19 

F20 

HIRXX(I) 


KMOM 


KPP 


KSIGN(I) 


NOMENCLATURE 

Balance components rotated and translated to the model 
(body) axis. 

Differential base pressure forces. 

Base force and moment tares. 

Final body axis components. 

Stability axis components. 

Wind axis components. 

Rotation from body axis to reference axis. 

Alternate reference axis coefficients. 

Reference axis coefficients. 

Base force and moment tare coefficients. 

Base pressure coefficients. 

Model (body) axis coefficients. 

Stability axis coefficients. 

Wind axis coefficients. 

Corrections for the effect of having a model with high 
restraints coupled with high interactions, where XX is the 
balance component (AF, SF, NF, RM, PM, YM) and I = 
balance number. 

Axial momentum correction term. 

= 0, no correction. 

= 1, applies nonblowing correction only and automatically 
computes FAMOM(I) 

= 2, applies nonblowing and blowing corrections 
A units conversion factor, initialized at 1. 

If PBASE is in PSF and PO is in PSI, KPP = 144.0 
If PBASE is in PSI and PO is in PSF, KPP = 0.00694 
If PBASE is differential (PBASE-PO), KPP = 0.0 
If PBASE is absolute, KPP = 1.0 (Standard). 

Constant for correcting balance quantities for grounding by 
wrong end, where I = balance number. 

KSIGN = 1 for normal balance attachment. 

KSIGN = -1 for grounding balance by wrong end. 
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SYMBOL 

D-7 

NOMENCLATURE 

k a,i 

COS(THETAO) * COS(PHIO) 

K A,2 

SIN(THETAO) 

K A,3 

COS(THETAO) * SIN(PHIO) 

L/D(I) 

Lift over drag ratio in the wind axis, where I = balance 
number. 

LS/DS(I) 

Lift over drag ratio in stability axis, where I = balance 
number. 

A 0 

Initial balance loads, roll moment, in. lb. 


Balance component quantity corrected for high interactions 
coupled with high model restraints, roll moment, in. lb. 

*4 

Balance component quantities corrected for balance 
orientation to gravity axis, roll moment, in. lb. 

METHOD 

Method to be used. 

MX(I) 

Final body axis rolling moment, in. lb., where I = balance 
number. 

MX(I,L) 

Balance rolling moment rotated (L = 1) and translated (L = 
2) to body axis, where I = balance number. 

MXREF'(I) 

Rolling moment rotated to reference axis, in. lb., where I = 
balance number. 

MXREF(I) 

Rolling moment translated to reference axis, in. lb., where I 
= balance number. 

MXS(I) 

Rolling moment in the stability axis, in. lb., where I = 
balance number. 

MXW(I) 

Rolling moment in the wind axis, in. lb., where I = balance 
number. 

MY(I) 

Final body axis pitching moment, in. lb., where I = balance 
number. 

MY(I,L) 

Balance pitching moment rotated (L = 1) and translated (L = 
2) to body axis, where I = balance number. 

MYREF'(I) 

Pitching moment rotated to reference axis, in. lb., where I = 
balance number. 

MYREF(I) 

Pitching moment translated to reference axis, in. lb., where 
I = balance number. > 
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SYMBOL 

MYS(I) 

MYW(I) 

MZ(I) 

MZ(I,L) 

MZREF'(I) 

MZREF(I) 

MZS(I) 

MZW(I) 


NF(I,J) 

NFO(I) 

NFT(I) 

NFTARE(I) 

NUBAL 

n o 

n 3 

n 4 


NOMENCLATURE 

Pitching moment in the stability axis, in. lb., where I = 
balance number. 

Pitching moment in the wind axis, in. lb., where I = balance 
number. 

Final body axis yawing moment, in. lb., where I = balance 
number. 

Balance yawing moment rotated (L = 1) and translated (L = 
2) to body axis, where I = balance number. 

Yawing moment rotated to reference axis, in. lb., where I = 
balance number. 

Yawing moment translated to reference axis, in. lb., where I 
= balance number. 

Yawing moment in the stability axis, in. lb., where I = 
balance number. 

Yawing moment in the wind axis, in. lb., where I = balance 
number. 

Initial balance loads, pitch moment, in. lb. 

Balance component quantity corrected for high interactions 
coupled with high model restraints, pitch moment, in. lb. 
Balance component quantities corrected for balance 
orientation to gravity axis, pitch moment, in. lb. 

Normal force, lbs., where I = balance number and J = 
correction number. 

Initial normal load, lbs., where I = balance number. 

Total normal load, lbs., where I = balance number. 

Normal weight tares, lbs., where I = balance number. 

Number of balances in the model. 

Initial balance loads, yaw moment, in. lb. 

Balance component quantity corrected for high interactions 
coupled with high model restraints, yaw moment, in. lb. 
Balance component quantities corrected for balance 
orientation to gravity axis, yaw moment, in. lb. 

Initial balance loads, normal force, lbs. 
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SYMBOL 



PBASE(II) 

PHI 

PHIB 

PHIB2 

PHIB3 

PHID 

PHID2 

PHID3 

PHIDX(I) 

PHIK 

PHIK2 

PHIK3 

PHIR 

PHI0,I 

PM 

PM(I,J) 

PMBASE(I) 

PMMAX 

PMO(I) 

PMT(I) 


NOMENCLATURE 

Balance component quantity corrected for high interactions 
coupled with high model restraints, normal force, lbs. 
Balance component quantities corrected for balance 
orientation to gravity axis, normal force, lbs. 

Base pressure, lbs/sq. in., where II = orifice number. 

Model Euler roll angle, degrees. 

Euler roll rotation angle between primary balance and 
model, degrees. 

Euler roll rotation angle between secondary balance and 
model, degrees. 

Euler roll rotation angle between tertiary balance and 
model, degrees. 

Roll deflection of primary balance, degrees. 

Roll deflection of secondary balance, degrees. 

Roll deflection of tertiary balance, degrees. 

Deflection roll angle constants, where X is balance 
component (A, S, N, R, P, Y) and I = balance number. 

Euler roll angle to account for knuckle and/or primary 
balance angles in relation to tunnel support, degrees. 

Euler roll angle to account for orientation of undeflected 
secondary balance in relation to primary balance, degrees. 
Euler roll angle to account for knuckle and/or tertiary 
balance angles in relation to tunnel support, degrees. 

Euler roll rotation angle between model (body) axis and 
reference axis, positive in same direction as PHIB, degrees. 
Wind off zero attitude of each balance, degrees, where I = 
balance number. 

Pitching moment, in. lb. 

Pitching moment, in. lb., where I = balance number and J = 
correction number. 

Base pitching moment, in. lb., where I = balance number. 
Maximum absolute value of pitch moment, in, lb. 

Initial pitching moment, in. lb., where I = balance number. 
Total pitching moment, in. lb., where I = balance number. 
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SYMBOL 

PMTAREI 

PSI 

PSIB 

PSIB2 

PSIB3 

PSID 

PSID2 

PSID3 

PSIDX(I) 

PSIK 

PSIK2 

PSIK3 

PSIR 

PSIU 

R(I,J) 

RGB 

RM 

RM(I,J) 

RMBASE(I) 

RMMAX 

RMO(I) 

RMT(I) 

RMTARE(I) 

SAREA(I) 


NOMENCLATURE 

Pitching weight tares, in. lb., where I = balance number. 
Model yaw angle, degrees. 

Euler yaw rotation angle between primary balance and 
model, degrees. 

Euler yaw rotation angle between secondary balance and 
model, degrees. 

Euler yaw rotation angle between tertiary balance and 
model, degrees. 

Yaw deflection of primary balance, degrees. 

Yaw deflection of secondary balance, degrees. 

Yaw deflection of tertiary balance, degrees. 

Deflection yaw angle constants, where X is the balance 
component (A,S,N,R,P,Y) and I = balance number. 

Euler yaw angle to account for knuckle and/or primary 
balance angles in relation to tunnel support, degrees. 

Euler yaw angle to account for orientation of undeflected 
secondary balance in relation to primary balance, degrees. 
Euler Yaw angle to account for knuckle and/or tertiary 
balance angles in relation to tunnel support, degrees. 

Euler yaw rotation angle between model (body) axis and 
reference axis, positive in same direction as PSIB, degrees. 
Tunnel sideflow angle, degrees. 

I'th row and J'th column in rotation matrix. 

Gravity to balance rotation matrix. 

Rolling moment, in. lb. 

Rolling moment, in. lb., where I = balance number and J = 
correction number. 

Base rolling moment, lbs., where I = balance number. 
Maximum absolute value of roll moment, in. lb. 

Initial rolling moment, in. lb., where I = balance number. 
Total rolling moment, in. lb., where I = balance number. 
Rolling weight tares, in. lb., where I = balance number. 
Reference area for balance coefficients. Normally wing 
area, sq. in., where I = balance number. 
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SYMBOL 

SF(I,J) 

SFO(I) 

SFT(I) 

SFTARE(I) 

TAREA 

TAREN 

TAREP 

TARER 

TARES 

TAREY 

THEDX(I) 

THETA 

THETAB 

THETAB2 

THETAB3 

THETAD 
THETA D2 
THETAD3 
THETAK 

THETAK2 

THETAK3 

THETAR 


THETAS 

THETAU 


NOMENCLATURE 

Side force, lbs., where I = balance number and J = correction 
number. 

Initial side load, lbs., where I = balance number. 

Total side load, lbs., where I = balance number. 

Side weight tares, lbs., where I = balance number. 

Axial momentum tare correction term. 

Normal momentum tare correction term. 

Pitching momentum tare correction term. 

Rolling momentum tare correction term. 

Side momentum tare correction term. 

Yawing momentum tare correction term. 

Deflection pitch angle constants, where X is the balance 
component (A,S,N,R,P,Y) and I = balance number. 

Model euler pitch angle, degrees. 

Euler pitch rotation angle between primary balance and 
model, degrees. 

Euler pitch rotation angle between secondary balance and 
model, degrees. 

Euler pitch rotation angle between tertiary balance and 
model, degrees. 

Pitch deflection of primary balance, degrees. 

Pitch deflection of secondary balance, degrees. 

Pitch deflection of tertiary balance, degrees. 

Euler pitch angle to account for knuckle and/or primary 
balance angles in relation to tunnel support, degrees. 

Euler pitch angle to account for orientation of undeflected 
secondary balance in relation to primary balance, degrees. 
Euler pitch angle to account for knuckle and/or tertiary 
balance angles in relation to tunnel support, degrees. 

Euler pitch rotation angle between model (body) axis and 
reference axis, positive in same direction as THETAB, 
degrees. 

Strut pitch angle, degrees. 

Tunnel upflow angle, degrees. 
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SYMBOL 

THETA0,(I) 

W 

x 

XBAR(I) 


XICH(I) 

XK 

XKCH(I) 

XREF 

y 

YBAR(I) 

YM 

YM(I,J) 

YMBASE(I) 

YMMAX 

YMO(I) 

YMT(I) 

YMTARE(I) 

YREF 


Z 

ZBAR(I) 

ZREF 


NOMENCLATURE 

Wind off zero attitude of each balance, degrees, where I = 
balance number. 

Weight tares. 

Distance of center of gravity to balance center, inches. 
Moment transfer distance measured in the body force axis 
system from the balance moment center to the desired 
moment center, positive in the direction of positive model 
thrust, side and normal force respectively, inches, where I = 
balance number. 

Intercept for momentum term, where I = balance number. 
Constants used in calculating momentum correction terms. 
Slope for momentum term, where I = balance number. 

Moment transfer distance. Measured relative to and in the 
same direction as XBAR, inches. 

Distance of center of gravity to balance center, inches. 

See XBAR. 

Yawing moment, in. lb. 

Yawing moment, in. lb., where I = balance number and J = 
correction number. 

Base yawing moment, lbs., where I = balance number. 
Maximum absolute value of yaw moment, in. lb. 

Initial yawing moment, in. lb., where I = balance number. 
Total yawing moment, in. lb., where I = balance number. 
Yawing weight tares, in. lb., where I = balance number. 
Moment transfer distance. Measured relative to and in the 
same convention as YBAR, inches. 

Initial balance loads, side force, lbs. 

Balance component quantity corrected for high interactions 
coupled with high model restraints, side force, lbs. 

Balance component quantities corrected for balance 
orientation to gravity axis, side force, lbs. 

Distance of center of gravity to balance center, inches. 

See XBAR. 

Moment transfer distance. Measured relative to and in the 
same convention as ZBAR, inches. 
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APPENDIX D 
Module D 

Balance Loads and Model Attitude 
A. Required Constants 

Required constants are defined in the nomenclatures. 

1. Primary balance deflection constants -A angle/A load 

PSIDA1 = A PSID/A AF(1,3) 

THEDA1 = A THETAD/A AF(1, 3) See related 

PHIDA1 = A PHID/A AF(1,3) items 2. and 3. 

PSIDN1 = A PSID/A NF(1,3) 

THEDN1 = A THETAD/A NF(1, 3) 

etc. 

2. Primary balance deflection angle names - PSID, THETA D, PHID. These 
names are optional as shown in item 3. However, they are suggested and 
extreme care should be used if changed since this is based on these names. 
No values are required for these angles since they are computed internally 
from the constants supplied under item 1 as follows: 


PSID = (PSIDA1)AF(1,3) + (PSIDN1)NF(1,3) 

+ (PSIDS1)SF(1,3) + (PSIDR1)RM(1,3) 

+ (PSIDP1)PM(1,3) + (PSIDY1)YM(1,3) (Eq. D-l) 

THETAD = (THEDA1)AF(1,3) + (Eq. D-2) 

PHID = (PHIDA1)AF(1,3) + (Eq. D-3) 


3. Input of items 1 and 2 - Deflection angle names and constants are input 
from C-card images (which may be modified) stored on magnetic storage 
disks. A maximum of six deflections is permitted. 


65 


Therefore, the six values assigned in the yaw plane (PSI) for example 
are PSIDA1, PSIDS1, PSIDN1, PSIDR1, PSIDP1, and PSIDY1 as 
defined in item 1. 


Input of rotations from gravity to primary balance - Rotations from 
gravity to primary balance axis system (see Figure D-l(a) to D-l(e)) are 
input from the R-card image names stored on magnetic storage disks. 

Secondary balance deflection constants - A angle/A load 


PSIDA2 

= A PSID2/A AF(2,3) 


THEDA2 

= A THETAD2/A AF(2,3) 

See related 

PHIDA2 

= A PHID2/A AF(2,3) 

Items 6. and 7. 

PSIDN2 

= A PSID2/A NF(2,3) 


THEDN2 

= A THETAD2/A NF(2,3) 


etc. 



Secondary balance deflection angle names - 

PSID2, PSID3, THETAD2, 

PHID2. These 

names are optional as shown in item 7. However, they are 

suggested and 

extreme care should be used if changed since this 


description is based on these names. No values are required for these 
angles since they are computed internally from the constants supplied 
under item 5 as follows: 

PSID2 = (PSIDA2)AF(2,3) + (PSIDN2)NF(2,3) 

+ (PSIDS2)SF(2,3) + (PSIDR2)RM(2,3) (Eq. D-4) 

+ (PSIDP2)PM(2,3) + (PSIDY2)YM(2,3) 

THETAD2 = (THEDA2)AF(2,3) + ( Eq . D _ 5) 

PHID2 = (PHIDA2)AF(2,3) + ( Eq . D _ 6) 

Input of items 5. and 6. - Deflection names and constants are input from 
C-card image names stored on magnetic disks. Six is the maximum 
number of deflections permitted. 
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8. Tertiary balance deflection angles are handled in a manner similar to 
primary and secondary balance constants. 

9. Input of rotations (THETAK2, PSIK2, THETA D2, etc.) from primary 
balance to secondary balance - Rotations from the primary balance to the 
secondary balance are input from R-card images stored on magnetic disks. 
See Figure D-l(f). 

10. Wind-off-zero attitude of each balance - Input PHIO, THETAO, from card 
images stored on magnetic disks for each balance. This option is normally 
used as a result of problems associated with option 2. It is also used when 
data zeros are not used in the force data reduction scheme. If data zeros 
are not taken and values are not input from the disk, PHIO = THETAO = 0 
is assumed. See Figure D-l(g). 

11. Weight tares and attitude loads - Tares are determined automatically 
from a 700 series weight-shift run made immediately before each model 
configuration tunnel run. Do not input W, X, Y, Z, W(AF), W(SF), . . . ., 
etc. 


12. HIRAFI, HIRNFI, HIRSFI, HIRPMI, HIRYMI, HIRRMI where I = balance 
number - These constants correct for the effect of having a model with 
high restraints (HIR) coupled with a balance with high interactions (AF, 
NF, etc.). Thus, the name HIRAFI, HIRNFI, etc. These constants are 
obtained for each balance component by the following equation. 


HIR xx(I) = 


Tunnel balance xx calibration 
xx span check 


(Eq. D-7) 


where xx = balance component 

Note that when this correction is applied, the balance spans should be 
used in the standard program for quantities (EU) and not in-tunnel 
calibration. These constants are input from the C-card images stored on 
the magnetic disks for each balance. 
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13. KSIGN(I) - Constant for correcting balance quantities for grounding by the 
wrong end, where I = balance number. As shown in Figure D-2, grounding 
the balance by the wrong end ("A" cases) rather than the taper end results 
in a change of each balance component sign. Therefore 

KSIGN(I) = 1 for normal balance attachment 
KSIGN(I) = -1 for grounding balance by wrong end. 

14. THETAU - Tunnel upflow angle, see Figure D-3. 

15. PSIU - Tunnel sideflow angle, see Figure D-3. 

16. Input of items 14. and 15. - THETAU and PSIU are the required rotations 
for the wind-to-gravity transformation and are input from the T-card 
images (tables as function of MACH) stored on magnetic disks. 

17. Euler yaw, pitch and roll rotation angles (PSIB, THETAB, PHIB) between 
balance and model, are shown in Figure D-4(a). 

18. Input of - PSIB(I), THETAB(I), and PHIB(I) - Required rotations for the 
balance-to-model transformation are input from C-card images stored on 
magnetic disks. 

19. XBAR(I), YBAR(I), ZBAR(I) - Moment transfer distances are measured in 
the body force axis system from the balance moment center to the 
desired moment center, positive in the direction of positive model thrust, 
side and normal force, respectively (see Figure D-4(b)). Input from C- 
card images stored on magnetic disks, where I = balance number. 

20. ARPB(II,K) - Areas or momentum arms * areas used with PBASE(II) for 
computing base force and base moment tares, where II = orifice number. 
Use care to insure proper tare force signs. Area and arm units must be 
consistent with units of base pressures and balance components. 
ARB(H,K) is the same but for the second balance. ARP(II,K) is the same 
but for the third balance. 
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21. Input of item 20. - Areas and arm x areas are input from C-card images 
stored on magnetic disks. A maximum of 20 may be used. 

22. KPP - Units conversion factor, initialized at 1. 

If PBASE is in PSF and PO is in PSI, KPP = 144 
If PBASE is in PSI and PO is in PSF, KPP = .0069444 
If PBASE is differential (PBASE-PO), KPP = 0 
If PBASE is absolute, KPP = 1 (standard) 

Input from C-card images stored on magnetic disks if not equal to 

1 . 0 . 

23. Input of items - PSIR, THETAR, and PHIR are the required rotations for 
the model (body) to reference axis transformation and are input from C- 
card images stored on magnetic disks. 

24. XREF, YREF, ZREF - Moment transfer distances are measured relative 
to and in the same convention as XBAR, YBAR and ZBAR. Input from C- 
card images stored on magnetic disks. 

B. Test for Balance Loads and Model Attitudes 

If NUBAL = 0, skip module D. 

C. Balance Component Naming System 

1. In general, the balance component naming system follows the format of 
WX(Y,Z), where 

WX = component name is as follows: 

AF = Axial force 
NF = Normal force 
SF = Side force 
PM = Pitching moment 
YM = Yawing moment 
RM = Rolling moment 
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Y = balance number associated with component 

1 = 1st balance 

2 = 2nd balance 
etc. 

Z - number of corrections applied to component 
(uncorrected quantity = 1). 


Uncorrected Balance Quantities 

1* Signs on component quantities are uncorrected and thus are a strict 
function of model-balance orientation and the manner in which the model- 
balance attachment is made. Figure D-2 provides sketches showing the 
eight most frequent cases of model-balance orientation and the 
corresponding component signs. Each case is shown for grounding the 
balance taper end and for grounding the balance opposite end ("A" cases). 

2. For normal NASA type balances, the component quantities are obtained 
directly from the standard program for quantities. The balance 
components for this type of balance are always named as follows: 

Axial force 
Normal force 
Side force 
Pitch moment 
Yaw moment 
Roll moment 

where I = balai 


- AF(I,1) 

- NF(I,1) 

- SF(I,1) 

- PM(I,1) 

- YM(I,1) 

- RM(I,1) 


C«] 


nee number 


(Eq. D-8) 


For TASK type balances, the component quantities are also obtained 
directly from the standard program for quantities (EU), but additional 
equations must be supplied since axial force and rolling moment are 
generally the only two components obtained directly with TASK type 
balances. The following equations and names are suggested for the 
engineering units program. The following equations assume the axes 
origin is at the center of the balance. 
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NF(I,1) = N1(I,1) + N2(I,1) ( Eq . D _g) 

PM(I,1) = N1(I,1) - N2(I,1) ( Eq . D _ 10) 

SF(I,1) = SI (1,1) + S2(I,1) ( Eq . D _ n) 

YM(I,1) = S1(I,1) - S2(I,1) (Eq. D-12) 


The names shown for the final quantities are mandatory. 

E. Tunnel Support Pitch Angle 

The tunnel support pitch angle is used in gravity to balance transformations. 

1. The tunnel support pitch angle is THETAS. See Figure D-l(a). 

2. THETAS is computed in the standard program for quantities. It may be 
obtained from the strut helipot or from a "dangle" meter in the model or 
input as a constant. 

F - Balance Quantities Corrected for Interactions. Weight Tares and Momentum 

Tares 

1. Balance component quantities corrected for interactions are named as 
follows: 


Axial force 

AF(I,2f 

Normal force - 

NF(I,2) 

Side force 

SF(I,2) 

Pitch moment - 

PM(I,2) 

Yaw moment - 

YM(I,2) 

Roll moment - 

_RM(1,2)_ J 



(Eq. D-13) 


2. Balance component quantities corrected for high interactions coupled 
with high interactions coupled with high model restraints are named as 
follows: 
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Axial force - ! 

| 

*AF(I,3)** 

Normal force - 

NF(I,3) 

Side force 

SF(I,3) 

Pitch moment - 

PM(I,3) 

Yaw moment - 

YM(I,3) 

Roll moment - 

JtM(I,3)_ 



(Eq. D-14) 


3. 


Balance component quantities corrected for balance orientation to gravity 
axis, attitude loads and weight tares are named as follows: 


Axial force 

"AF(I,4)~ 

Normal force - 

NF(I,4) 

Side force 

SF(I,4) 

Pitch moment - 

PM(I,4) 

Yaw moment - 

YM(I,4) 

Roll moment - 

LRM(I,4)_ 



(Eq. D-15) 


4. Initial balance loads or weight tares are named as follows: where I = 

balance number. 


'AFO(I), NFO(I), SFO(I) 
JPMO(I), YMO(I), RM0(I)_ 



(Eq. D-16) 


5. Total balance loads (AF(I,1) + AFO(I), NF(I,1) + NFO(I), etc. 
are named as follows: 


AFT(I), NFT(I), SFT(I) ^ 
PMT(I), YMT(I), RMT(Ii. 



(Eq. D-17) 


6. First order interactions £ire represented by a matrix Cl; second order 
interactions are represented by a matrix C2. 
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7. Attitude weight tares are named as follows: 


AFTARE(I), NFTARE(I), SFTARE(I) 
PMTARE(I), YMTARE(I), RMTARE(I) 

8. Constants required from the project engineer are: 


[ f tare] 


(Eq. D-18) 


a. For gravity-to-primary-balance rotations, see Figure D-l(e). 
For gravity-to-tunnel-strut rotation, see Figure D-l(a). 

THETAS is already supplied from E. 


For tunnel strut-to-undeflected-primary balance rotations, see 
Figure D-l(b) and D-l(c). 

PSIK, THETAK, PHIK 


For undeflected balance-to-deflected-balance rotations, see Figure 
D-l(d). 

PSIDA1, THEDA1, PHIDA1 
PSIDS1, THEDS1, PHIDS1 
PSIDN1, THEDN1, PHIDN1 
PSIDR1, THEDR1, PHIDR1 
PSIDP1, THEDP1, PHIDP1 
PSIDY1, THEDY1, PHIDY1 

b. Primary- to-secondary-balance rotations 

For primary balance-to-undeflected-secondary balance rotations, 
see Figure D-l(f). 

PSIK2, THETAK2, PHIK2 
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Undeflected secondary balance-to-deflected-secondary balance 
rotations (with respect to primary balance). 

PSIDA2, THEDA2, PHIDA2 
PSIDS2, THEDS2, PHIDS2 
PSIDN2, THEDN2, PHIDN2 
PSIDR2, THEDR2, PHIDR2 
PSIDP2, THEDP2, PHIDP2 
PSIDY2, THEDY2, PHIDY2 

The third balance is similar to the above but with the number 3 
replacing the number 2 in the second balance. 

For wind-off-zero attitude of each balance (See Figure D-l(a)) 

PHIO, I, THETAO, I, 

c. High restraint and interaction constants 

HIRAFI, HIRNFI, HIRSFI 
HIRPMI, HIRYMI, HIRRMI 


The following description on correcting balance quantities for interactions 
and weight tares does not provide the exact equations for computing 
corrected balance quantities. The PAB balance check point program or 
the contractor's user manual must be consulted for these. However, this 
does provide the general outline for computing corrected balance 
quantities. 


Determine uncorrected total loads, (fut] 


[fut] = (fi)+ (fcT) = 


~AF(I,1) + AF0(I)~ 
SF(I,1) + SFO(I) 
NF(I,1) + NFO(I) 
RM(I,1) + RMO(I) 
PM(I,1) + PMO(I) 
_YM(I,1) + YM0(I)_ 


(Eq. D-19) 


Correct for interactions 
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a. [FUf] = (CQ* [FT] + 0 2 D* (Eq. D-20) 


where 



are balance interaction constants 


b. Therefore 


[FT] = [C] _1 *(FUT] - [Cj]" 1 *(p0* 0P] (Eq. D-21) 


Compute corrected delta balance loads, 




'AF(I,2r 


^AFTa) - AFoar 

SF(I, 2) 


SFT(I) - SFO(I) 

NF(I,2) 

= 

NFT(I) - NFO(I) 

RM(I,2) 


RMT(I) - RMO(I) 

PM(I,2) 


PMT(I) - PMO(I) 

JM(I,2£ 


JfMT(I) - YM0(I)_ 


(Eq. D-22) 


10. Correct forces and moments for high model restraints coupled with high 
balance interactions 



~AF(I,3) - 


'AF(I,2) + (HIRAF)AF(I,1) " 

SF(I,3) 


SF(I,2) + (HIRSF)SF(I,1) 

NF(I,3) 

= 

NF(I,2) + (HIRNF)NF(I,1) 

RM(I,3) 


RM(I,2) + (HIRRM)RM(I,1) 

PM(I,3) 


PM(I,2) + (HIRPM)PM(I,1) 

_YM(I,3^ 


_YM(I,2) + (HIRYM)YM(I,1) 


(Eq. D-23) 


11. Depending on the value of the constant KMOM, balance components are 
further corrected for balance /b ellows interactions and momentum flow 
effects. 
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If KMOM = 0, (Eq. D-24) 

no further balance corrections are applied and equations D-29 to D- 
35 are skipped. 

If KMOM = 1, 

nonblowing balance corrections are applied and FAMOM(I) is 
atuomatically computed. 

APCH = 0.0 
AMOM(I) = 0.0 

Equations D-28 to D-38 are executed. 

FJCON/FI = f(PTJ/PO) Table lookup 

FAMOM(I) = AF(I,4) - FI [FJCON/FI] (Eq. D-25) 

The values of FJCON/FI are obtained from an input table which results 
from averaged Stratford choke nozzle data obtained over many years. 
Typical table values are given below: 


PTJ/PO 

FJCON/FI 

PTJ/PO 

FJCON/FI 

1.0 

0.0 

5.0 

0.9700 

1.3 

0.9820 

6.0 

.9600 

1.5 

.9905 

7.0 

.9500 

2.0 

.9960 

8.0 

.9425 

3.0 

.9920 

10.0 

.9300 

4.0 

.9815 

14.0 

.9125 

4.5 

.9760 




A maximum of 15 values can be input to the computer as a T table. 

If KMOM = 2 and PCH(I) < 25, (Eq. D-26) 

the jet is assumed to be off and only nonblowing balance corrections 
are applied. 


APCH = 0.0 
AMOM(I) = 0.0 

Equations D-28 to D-38 are executed. 


D-25 


If KMOM = 2 and PCH(I) >25, 

the jet is assumed to be operating and all balance corrections are 
applied. 

APCH = PCH(I) 

Equations D-27 to D-38 are executed. 

Double second-order curve capability for computation of AMOM(I). I = 
balance number 

If APCH < XK ?3 j 

AMOM(I) = XKCH(I) * APCH + XICH(I) + XK„. * APCH 2 

74,1 

If APCH>XK then 

AMOM(I) = XKCH(I + 3) * APCH + XICH(I + 3) + XK * APCH 2 

75,1 

(Eq. D-27) 

"AF(I,3) 

SF(I,3) 

NF(I,3) 

RM(I,3) 

PM(I,3) 

>YM(I,3)J 

Balance/bellows interactions and momentum flow effects on the balance 
are computed after high restraint corrections. 
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TAREA 


TAREN 


TAREP 


TARES 


TAREY 


TARER 


AMOM(I) + XKj j * SF + XK 2 l * FN + XK 3 j * RM (Eq. D-29) 
+ XK A l * PM + XK 5 ^ * YM + XKg j * APCH 


(APCH - XK g r ) * [XKg j * SF + XK J0 j * FN 

* XK 11,I * R 1 * 1 + XK 12,I* PM + XK 13,I * YM ] 

* (APCH - XK U I ) • XK 15j , 4 XK 16>i « SF 

+ XK 17,I * RM + XK 18,I * PM + XK X9,I * ™ 


(APCH - XK,. .) * fxK . » SF * XK„ . * FN 

* XK 24,I * RM ’ + XK 25,I * PM + XK 26,I * ™3 

+ (APCH - XK„ .) « XK-. . + XK,. . * SF 
tXK 30,I ,FN+,XK 31,l*R MtXK 32,I* YM 


(A p CH - XK 34 I ) • [xk 35>i * SF 4 XK 36j , • FN 

* XK 37,I * RM + XK 38,I * PM + XK 39,I * YM ] 

* (APCH - XK 40 i ) • XK 41> , 4 XK,, , • FN 

* XIi 43,I * RM + XK 44,I * PM + XK 45,I * YM 


(APCH - XK ) * [XK * SF t XK JS ♦ FN 
7 XK 50i • RM4 XK 51j] • PM 4 XK^, -YM] 

* (APCH - XK 53>i ) • XK 54>i * XK 55j , • SF 

+ XK_ C T * FN + XK _ . * RM + XK co . * PM 

objl 5 1 jl 5o 

+ XK 59,I 

<APCH * XK 60,I> * C XK 61,I * SF + XK 62,I * FN 

* XK 63; i • RM 4 XK 64>i • PM 4 XK 65j , * YM] 

4 (APCH - XK, 6iI ) • XK 6?>i 4 xk 68>i • SF 

* XK 69,I * ™ + XK 79,I * PM * Xli 71,I * YM 

* XIi 72,I 


(Eq. D-30) 


(Eq. D-31) 


(Eq. D-32) 


(Eq. D-33) 


(Eq. D-34) 
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"AF(I,3r 


"AF(I,3) - TAREA" 

SF(I,3) 


SF(I,3) - TARES 

NA(I,3) 


NF(I,3) - TAREN 

RM(I,3) 


RM(I,3) - TARER 

PM(I,3) 


PM(I,3) - TAREP 

_YM(I,3)_ 

k 

_YM(I,3) - TAREY_ 


12. Perform gravity-to-balance transformations. 

Let (R0 denote specific Euler transformation matrixes 

[ Fba D_ PstrutD Pknuckle] PdeflectionsD (?g 

" L r gbJ (Jg] 

where 

C F baD~ vector representing balance quantities 
?g] = vector representing balance quantities 
G*Gb) = gravity-to-balance axis transformation 
13. Determine weight tares (attitude loads) 


"a FT A RE " 


"W(sin 0 - sin 0J "t 

SFTARE 

= 

o ^ 

W(cos 0^ sin - cos 0^ sin <j>g) 

NFTARE 


-W(cos 0^ cos <j>^ - cos 0q COS 4»g) 

RMTARE 


SFTARE(Z) - NFTARE(Y) 

PMTARE 

I 

; 

AFTARE(Z) + NFTARE(X) 

JYMTARE_ 


.SFTARE(X) + AFTARE(Y) 


Correct for weight tares (attitude loads) 


(Eq. D-35) 


J (Eq. D-36) 

in balance axis, 
in gravity axis, 
matrix. 

(Eq. D-37) 
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(F43 = (F3] - (FTARE^ = 


'AF(I,4F 


"AF(I,3) - AFTARE(I)" 

SF(I,4) 


SF(I,3) - SFTARE(I) 

NF(I,4) 

= 

NF(I,3) - NFTARE(I) 

RM(I,4) 


RM(I,3) - RMTARE(I) 

PM(I,4) 


PM(I,3) - PMTARE(I) 

_YM(I,4) 


_YM(I,3) - YMTARE(Ii. 


(Eq. D-38) 


G. Balance Quantities Corrected for Method of Attachment 

1. Balance component quantities corrected for method of attachment are 
named as follows: 


F5I = 


AF(I,5) 

SF(I,5) 

NF(I,5) 

RM(I,5) 

PM(I,5) 

YM(I,5)J 


(Eq. D-39) 


Where I = balance number. 


2. The constant required from the project engineer is KSIGN(I). 


(j?5) = KSIGN *(F4) = 


*F(I,5f 


SF(I,5) 


NF(I,5) 

= 

RM(I,5) 


PM(I,5) 


YM(I,5i 



’KSIGN(I) * AF(I,4)" 
KSIGN(I) * SF(I,4) 
KSIGN(I) * NF(I,4) 
KSIGN(I) * RM(I,4) 
KSIGN(I) * PM(I,4) 
KSIGN(I) * YM(I,4)_ 


(Eq. D-40) 


H. Angle of Attack and Sideslip Angle 

1. The following definitions denote various transformation matrixes which 
are obtained from given orders of Euler rotation angles. 
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wind-axes-to-gravity-axes transformation matrix 
gravity-axes-to-balance-axes transformation matrix. This 
matrix is established from rotation angles supplied in section 
F, therefore 

R strut R knuckle Reflection ( E< 1* D-41 ) 

balance axes-to-model axes transformation matrix 

2. The constants required from the project engineer are THETA U, PSIU, 
PSIBI, THETABI AND PHIBI. 

For wind-to-gravity rotation angles, see Figure D-3. 

For balance-to-model rotation angles, see Figure D-4(a). 

The matrix > which transforms a vector in the wind axis system to the model axis 

system, may now be computed by a yaw, pitch, and roll rotation. The result is the final 
rotation matrix from the wind axes to model axes. 



ti 

C R B10C E GeDC R W<l) * 

" W ll w 12 W 13 _1 

W 21 w 22 w 23 



w W 31 w 32 w 33 - 

ii 

CVD [VDCV'f 

) 


(Eq. D-42) 


(Eq. D-43) 
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C r wm] = 



'1 

0 O' 


'cos 0 0 -sin 0" 

= 

0 

cos 4> -sin <t> 


0 10 


_0 

sin <t> cos <t>_ 


_sin 0 0 cos9_ 

= 


” cos 0 

0 

-sin 6 



-sin <j) sin 0 cos 4> 

-sin 4> cos 6 



^cos <t> sin 0 sin <t> 

cos <J> cos 0_ 


cos 0 cos 


-sin ^ cos 0 


cos 'J' 
sin ^ 
. 0 


cos 'P 
sin 'I' 
L 0 


-sin 'J' 0 ' 

cos ^ 0 

0 1 . 

-sin 4 1 0 ' 

cos ^ 0 

0 1 . 


-sin <|» sin 0 cos ^ + cos 4> sin \|> 
_ cos (J> sin 0 cos + sin 4> sin 


sin 4> sin 0 sin + cos <l> cos 
-cos <|> sin 0 sin ^ + sin <l> cos ^ 


-sin 9 

-sin <t> cos 9 
cos <J> cos 9 

(Eq. D-44) 


0 Pitch angle 

<t> Roll angle 

i|> Yaw angle 


Using the definitions shown in Figure D-5 and the above information 


ALPHA = TAN -1 


W 31 
(rr 1 -) 


w 


11 


Note that for 4>= 0°, a = 0 

PSI = SIN -1 (w 21 ) 

BETA = -PSI 

THETA = SIN _1 (-W 13 ) 
-1 w 23 

PHI = TAN 1 (- — ^-) 
w 

W 33 


(Eq. D-45) 


(Eq. D-46) 
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I. 

1. Balance components rotated to the model (body) axis are named as 
follows: 

Axial - FA(I,1) 

Side - FY(I,1) 

Normal - FN(I,1) 

Roll - MX(I,1) 

Pitch - MY(I,1) 

Yaw - MZ(I,1) 

2. Balance components rotated and translated to the model (body) axis are 
named as follows: 

Axial - FA(I,2) 

Side - FY(I,2) 

Normal - FN(I,2) 

Roll - MX(I,2) 

Pitch - MY(I,2) 

Yaw - MZ(I,2) 

3. The constants required from the project engineer are XBAR, YBAR and 
ZBAR. (See Figure D-4.(b)) 

The matrix (j* BI yQis used to transform the components in the balance axis to 
the model (body) axis system as follows: 


(Eq. D-47) 


(Eq. D-48) 
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or 


FA(i,ir 


^b n AF(I,5) + b 12 SF(I,5) + b 13 NF(I,5) -j 

FY(I,1) 


b 21 AF(I,5) + b 22 SF(I,5) + b 23 NF(I,5) 

FN(I,1) 


b 31 AF(I,5) + b 32 SF(I,5) + b 33 NF(I,5) 

MX(I,1) 


b 11 RM(I l 5) - b 12 PM(I,5) + b 13 YM(I,5) 

MY(I,1) 


-b 21 RM(I,5) + b 22 PM(I,5) - b 23 YM(I,5) 

MZ(I,1)_ 


^ b 31 RM(I,5) - b 32 PM(I,5) + b 33 YM(I,5) 


(Eq. D-49) 


The components are then translated as follows 


FA(I,2) 

FY(I,2) 

FN(I,2) 

MX(I,2) 

MY(I,2) 

w MZ(I,2) 


FA(I,1) 

F Y(I,1) 

FN(I,1) 

MX(I,1) + FN(I,1) * YBAR - FY(I,1) * ZBAR 
MY(I,1) - FN(I,1) * XBAR - FA(I,1) * ZBAR 
MZ(I,1) - FY(I,1) * XBAR - FA(I,1) * YBAR_ 


(Eq. D-50) 


J. Pressure Corrections to Body Axis Components 

1. Base and/or cavity pressures are obtained from the standard program for 
quantities and are named PBASE(II). Where II = orifice number. 


2. Tunnel static pressure is computed in module A and is named PO. 

3. Base force and moment tares are named as follows: 


Axial - FABASE(I) 
Side - FYBASE(I) 
Normal - FNBASE(I) 
Roll - RMBASE(I) 
Pitch - PMBASE(I) 
Yaw - YMBASE(I) 
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4. Final body axis components, corrected for base tares, are named as 
follows: 

Axial - FA(I) 

Side - F Y(I) 

Normal - FN(I) 

Roll - MX(I) 

Pitch - MY(I) 

Yaw - MZ(I) 

Note that axial force is not corrected for internal (duct) axial force. 

5. The constants required from the project engineer are ARPB(II,K) and 
KPP. 

To determine differential base and cavity pressures 

A PBASE(II) = PBASE(H) - (fro * (KPP))), (Eq. D-51) 

Noting that a positive differential pressure acting on the base of a model causes 
a thrust, then base pressure force and moment tares are defined as follows: 

FABASE(I) = -Z (APBASE(II))*(ARPB(II,1)) (Eq. D-52) 

11=1 

FYBASE(I) = l [A PBASE(II)) * [aRPB(II, 2)) (Eq. D-53) 

11=1 


n r 

FNBASE(I) = Z [a PBASE(II)j * (aRPB(II, 3)) (Eq. D-54) 

11=1 
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RMBASE(I) = 2 [APBASE(II))*(ARPB(II,4)]) 
11=1 

PMBASE(I) = 2 |APBASE(II)j * QvRPB(II,5)J 
11=1 

YMBASE(I) = 2 (APBASE(Il[|*[ARPB(n,6g 
11=1 


^FA(I)~ 


"FAd^r 


'fabase(i)^ 

F Y(I) 


FY(I,2) 


FYBASEd) 

FN(lj 

= 

FN(I,2) 

- 

FNBASE(I) 

MX(I) 


MX(I,2) 


RMBASE(I) 

M Y(I) 


MY(I,2) 

[ 

PMBASE(I) 

^MZd)^ 


JMZ(I,2)_^ 


jrMBASE(I)_ 


Stability Axis Components 

1. Force and moment components in the stability axis are called 


Drag - 

FDS(ir 

Side 

FYS(I) 

Lift 

FLS(I) 

Roll - 

MXS(I) 

Pitch - 

MYS(I) 

^Yaw - 

MZS(I)_j 


where I = balance number. 


Note that drag is not corrected for internal (duct) drag. 


fds(i) = [fa(i])* (cos(alpha)) + [fn(i^* (sin(alpha)) 

FYS(I) = FY(I) 

FLS(I) = [FN(l]|* (C0S(ALPHA3 -[fA(I)) ‘(SIN(ALPHA)) 


(Eq. D-55) 


(Eq. D-56) 


(Eq. D-57) 


(Eq. D-58) 


(Eq. D-59) 
(Eq. D-60) 
(Eq. D-61) 
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MXS(I) = (MX(I)] * (COS(ALPHA)) + (MZ(l[| *(SIN(ALPHA[) (Eq. D-62) 

MYS(I) = MY(I) (Eq.D-63) 

MZS(I) = (MZ(I)3*(COS(ALPHA3-[MX(I)3*0IN(ALPHA3 (Eq. D-64) 


L. Wind Axis Components 

1. Force and moment components in the wind axis are called 


Drag - FD(I) 

Crosswind - FC(I) 

Lift - FL(I) 

RoU - MXW(I) 

Pitch - MYW(I) 

Yaw - MZW(I) 


Note that drag is not correct for internal (duct) drag. 


FD(I) = (F DS(lj] * (COS(BETA)] - (f YS(I)] *(siN(BETA)] (Eq. D-65) 

FC(I) = (FYS(I)) * (cOS(BETA)j + (FDS(lf) *(SIN(BETA)) (Eq. D-66) 

FL(I) = FLS(I) (Eq. d_ 67) 

MXW(I) = (mXS(I)) * (COS(BETA)) + (mYS(I)1 *fsiN(BETA)] (Eq. D-68) 

MYW(I) = (MYS(I)] * |oS(BETA)] - (mXS(I)J *(siN(BETA0 (Eq. D-69) 

MZW(I) = MZS(I) ( E q. D _ 70) 


M. Alternate Reference Axis Components 

1. Body axis components rotated and translated to an arbitrary reference 
axis system are called 

Axial - FAREF(I) 

Side - FYREF(I) 

Normal - FNREF(I) 

Roll - MXREF(I) 

Pitch - MYREF(I) 

Yaw - MZREF(I) 
where I = balance number. 
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Note that axial force is corrected for internal (duct) axial force. 


2. The transformation matrix for model axis to reference axis rotations is 
defined as ^R M ^j. 

3. The constants required from the project engineer are PSIR, THETAR, 
PHIR, XREF, YREF, ZREF and SAREAI where I = balance number for 
model-{body)-to-reference axis rotations. 

4. CAI is from module E. 


The matrix (r mr J is used to transform the components in the model (body) axis 
to a reference axis system as follows: 


FA(I)' = FA(I) - CAI * QO * SAREA(I) 


(Eq. D-71) 


FAREF(I)' 

F YREF(I)' 
FNREF(I)' 

-MXREF(I)' 

MYREF(I)' 

-MZREF(I)' 


FAREF(I)' 

FYREF(I)' 

FNREF(I)' 

MXREF(I)' 

MYREF(I)' 

MZREF(I)' 


[ r mr] 



FA(I)' 

FY(I) 

FN(I) 

-MX(I) 

MY(I) 

-MZ(I) 


"m n FA(IV + m 12 FY(I) + m 13 FN(I) “ 
m 21 FA(I) + m 22 FY(I) + m 23 FN(I) 
m 31 FA(I)’ + m 32 FY(I) + m 33 FN(I) 
m n MX(I) - m 12 MY(I) + m 13 MZ(I) 
-m 21 MX(I) + m 22 MY(I) - m 23 MZ(I) 
.m 31 MX(I) - m 32 MY(I) + mggMZd)^ 


(Eq. D-72) 


(Eq. D-73) 


(Eq. D-74) 


The components are now translated as follows 
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"FAREF(lf 


Tfaref(i)' 

FYREF(I) 


FYREF(I)' 

FNREF(I) 

_ 

FNREF(I)' 

MXREF(I) 


MXREF(I)' + FNREF(I)' * YREF - FYREF(I)' * ZREF 

MYREF(I) 


MYREF(I)' - FNREF(I)' * XREF - FAREF(I)' * ZREF 

JVlZREFd^ 


^MZREF(I)' - FYREF(I)’ * XREF - FAREF(I)' * YREF 


(Eq. D-75) 

N. Base Force and Moment Tare Coefficients 

1. Base force and moment tare coefficients are called 

"Axial - CABASE(I) " 

Side - CYBASE(I) 

Normal - CNBASE(I) 

RoU - CRMBASE(I) 

Pitch - CPMBASE(I) 

Jaw - CYMBASE(I)^ 

where I = balance number. 

2. Free-stream dynamic pressure is defined in module A and is called QO. 

3. The constants required from the project engineer are SAREA(I), 
CHORD(I), and BSPAN(I). 


CABASE(I) ' 


~FABASE(I) 

CYBASE(I) 


FYBASE(I) 

CNBASE(I) 

1 

FNBASE(I) 

CRMBASE(I) 

j^QO * SAREAdrj 

[RMBASEd)/BSPANdfi 

CPMBASE(I) 


[PMBASEd)/CHORDd)j 1 

_CYMBASE(IL 


w &MBASE(I)/BSPAN(I)jJ 


(Eq. D-76) 
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O. Base Pressure Coefficients 

1. Base pressure coefficients are called CPBASEdI) 

CPBASEdI) = ^[A PBASEdlj] (Eq. D-77) 

where II = orifice number. 

P. Model (Body) Axis Coefficients 

1. Model (body) axis coefficients are called 

Axial - CAd) 

Side - CY(I) 

Normal - CNd) 

Roll - CMXd) 

Pitch - CM Y(I) 

Yaw - CMZ(I) 

where I = balance number. 

2. CAI is from module E. 


CAd) 


FAd) 


CAI 

CY(I) 


FY(I) 


0 

CN(I) 

1 

FN(I) 


0 

CMXd) 

" p<3]*|3AREA(l}] 

MX(I)/BSPAN(I)] 

- 

0 

CM Yd) 


|MYd)/CHORD(I)] 


0 

CMZd) 

* 

faZ(I)/BSPAN(I)] 


0 


(Eq. D-78) 
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Q. Stability Axis Coefficients 

1. Stability axis coefficients are called 


Drag - CDS(I) 

Side - CYS(I) 

Lift - CLS(I) 

Roll - CMXSG) 

Pitch - CMYS(I) 

Yaw - CMZSd) 


where I = balance number. 


2. CDIS is from module E. 


CDSd) 

CYSd) 

CLSd) 

CMXSG) 

CMYSG) 

CMZSG) 


FDSG) 

FYSft) 

FLSfl) 

[QCT]*[$AREAGJj jMXS(I)/BSPANGf| 

jMYSflVCHORDG^I 
|MZSd)/BSPANd]| 


R. Wind Axis Coefficients 

1. Wind axis coefficients are named 


Drag - CDG) 

Crosswind - CCG) 

Lift - CL ft) 

Roll - CMXWft) 

Pitch - CMYWft) 

Yaw - CMZWft) 


where I = balance number. 



2. CDI is from module E. 


— — 
CDd) 


FDd) 


CDI 

ccd) 


FC(I) 


0 

CLd) 

1 

FLd) 


0 

CMXW(I) 

“ p6]*[BAREAdl] 

jMXWd)/BSPANdj] 

- 

0 

CMYW(I) 


|MYW(I)/CHORD(lJ 


0 

CMZWd) 


[MZW(I)/BSPANd)] 


0 

!_ _ 


Alternate Reference Axis Coefficients 


(Eq. D-80) 


1. Reference axis coefficients are named 


Axial - CAREF(I) 
Side - CYREFO) 
Normal - CNREF(I) 
Roll - CMXREF(I) 
Pitch - CMYREFd) 
Yaw - CMZREFd) 


where I = balance number. 


CAREFd) 


FAREF(I) 

CYREF(I) 


FYREFd) 

CNREF(I) 

1 

FNREFd) 

CMXREFd) 

" £QCFJ*£S A R E A d TJ 

_MXREF(I)/BSPAN(I)] 

CMYREFd) 


MYREF(I)/CHORD(I)] 

CMZREFd) 


[MZREFd)/BSPA^(l)] 


(Eq. D-81) 

Miscellaneous Equations 

1. Base drag coefficient is called CDBASE(I). Where I = balance number. 


[c DBASE (I) = |cABASEd)|*|C0S(ALPHA)J+ JCNBASEd)j*[siN(ALPHA)J 


(Eq. D-82) 
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2. Lift-over-drag ratio in the stability axis is called LS/DSfl). 

LS/DSd) = CLSOVCDSd) (Eq. D-83) 

3. Lift-over-drag ratio in tbe wind axis is called L/D(I). 

L/D(I) = CL(I)/CD(I) (Eq. D-84) 

4. Lift coefficient squared is called CLSQRd). 

CLSQR(I) = [cLS(l)]* [CLS(I)] (Eq. D-85) 

U. Calculation of Initial Weight Tares and Attitude Load Constants 

1. The initial weight tares and attitude load constants may be obtained by 
either of two methods for each strain gage balance. 

a. Method I - Data obtained at an arbitrary series of pitch angles (2 £ 
number of pitch angles £ 30). This method cannot be used with a 
balance without an axial force component. 

b. Method II - Data obtained at an arbitrary series of roll angles (4 £ 

number of roll angles £ 30). Normally, the roll angles will be 0°, 

90°, 180°, and 270°. The roll angle must be specified in a digital 

channel with name PHIK. (Note that this method must be used for 
balances without an axial force component). This method cannot be 
used with a balance that does not have a rolling moment coefficient. 

V. Calculation of Initial Weight Tares and Attitude Load Tares (Strain Gage 
Balance) 

1. Calculate 

a * K A,1 = cos ** cos Bo (Eq. D-86) 

b. K A>2 = sin 0 o (Eq. D -87) 
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C * K A,3 = sin<t> ° cos 80 (Eq. D -88) 

2. Determine from balance deck number of components and what these 
components are. 

3. Determine maximum value of each equipment over entire tare run. 


a. 

FNMAXG) = ABS(NF(I,l))max 

(Eq. D-89) 

b. 

FAMAXd) = ABS(AF(I,l))max 

(Eq. D-90) 

c. 

FYMAXd) = ABS(SF(I,l))max 

(Eq. D-91) 

d. 

PMMAXd) = ABS(PMd,l))max 

(Eq. D-92) 

e. 

RMMAXd) = ABS(RMd,l))max 

(Eq. D-93) 

f. 

YMMAXd) = ABS(YM(I,l))max 

(Eq. D-94) 


4. Initialize initial weight tares and attitude load constants. 

a. Set A A= A N= A Y=0 (Eq. D-95) 

A m 1 = A m 2 = A n I = A n 2 = A *' 1 = Ajl 2 = 0 
x = y = z = 0 

b. Assume NFOG) = AFOft) = PMOG) = RMOtt) = YMOG) = SF0(I) = 0. 

(Eq. D-96) 

5. For each data point correct balance quantities for interactions. 

Determine uncorrected total loads, JfUt] 

AF (1,1) + AF0(I) 

SF (1,1) + SF0(I) 

NF0,1) + NFOG) (Same as Eq. D-19) 

RMG,1) + RMOG) 

PMa.i) + PMoa) 

YMd,l) + YMOd) 


[fut] . [F,> [fo] . 
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Correct for interactions 


a. [FUT3 = [CJ*(FT3 + [C 2 ]* C FP D (Same as Eq. D-20) 
and pjare balance interaction constants 

b. Therefore 


where (c/j 


[ft] = [Cj J 1 •[put] - [CJ 1 *[C0* [fp] 

(Same as Eq. D-21) 


Compute corrected delta balance loads, (j2^ 



~AF(I,2) " 


~AFT(I) - AFO(I) ~ 

SF(I, 2) 


SFT(I) - SFO(I) 

NF(I,2) 

= 

NFT(I) - NFO(I) 

RM(I,2) 


RMT(I) - RMO(I) 

PM(I,2) 


PMT(I) - PMO(I) 

^YM(I,2)_ 


^YMTd) - YMOtt)^ 


(Same as Eq. D-22) 


Correct forces and moments for high model restraints coupled with high 
balance interactions 



”AF(I,3)~ 


"AF(I,2) + (HIRAF)AF(I,1) ~ 

SF(I,3) 


SF(I,2) + (HIRSF)SF(I,1) 

NF(I,3) 

= 

NF(I,2) + (HIRNF)NF(I,1) 

RM(I,3) 


RM(I,2) + (HIRRM)RM(I,1) 

PM(I,3) 


PM(I,2) + (HIRPM)PM(I,1) 

YM(I,3) 


JfM(I,2) + (HIRYM)YM(I,1)^ 


(Same as Eq. D-23) 
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R GB " 


w. 


6. Determine balance rotation from gravity axis. 

a. Determine rotation matrix for each matrix. See first part of this 
module. 


b. Determine = product of each individual rotation 

c. Then: 


cos 0cos \p 

-sin <p sin 0 cos \p + cos <p sin ip 
. cos <j> sin 0 cos ip + sin <J> sin \p 


-sin ip cos 0 

sin 4> sin 0 sin ip + cos <p cos ip 
-cos <p sin 0sin \p + cps (pcos ip 


-sin 0 - 

-sin <p cos 0 
cos <pcos 0_ 


R(l,l) 

R(l,2) 

R(l,3) 

R(2,l) 

R(2,2) 

R(2,3) 

. R(3,l) 

R(3,2) 

R(3,3) 


(Eq. D-97) 
(Eq. D-98) 


e. Calculate 

PSI = SIN _1 (w 21 ) 
BETA = -PSI 


THETA = SIN" x (-W 13 ) 


PHI = TAN -1 ( 


w 


w 


23_ 

33 


(Same as Eq. D-46) 


Calculation of Attitude Load Constants by Method I 

1. Solve following matrix equation using a least squares technique (MINFIT 
routine) for A A. 
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(-R(1,3)-K a 2 )j 
(-R(1,3)-K a ; 2 ) 2 


AA 


(-R(l,3)-K A>2 ) k | 

where k is the number of data points .1 30 


(A,), 

« 3 ) 2 




(Eq. D-99) 


2. AN = AA + AW 


(Eq. D-100) 


AY = AN 


where AW is obtained from balance interaction deck 

3. If PMMAXfl) >YMMAX0) and 5RMMAX0) 

a. Solve following matrix equation using least squares technique for 
A m i and Anr» 2 . 


<K a j - R(3,3», 

Vl - R(3,3)) 2 
• 

(-R(1,3)-K a 2 ), 
(-R(1,3)-K a ; 2 ) 2 
• 


Am j 


(m 3>! 

(m 3 , 2 

• 

(K A> 1 - R(3,3)) k 

• 

<-R(l,3)- K A 2 ) k 


Am 2 


( m 3 ) k 


Am j 
“ AN 

(Eq. D-101) 

(Eq. D-102) 

Am- 


= AA - 

(Eq. D-103) 
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d. A = Am 2 (Eq. D-104) 

e. An^ = Anij (Eq. D-105) 

f. If YMMAX(I) > RMMAX(I) solve the following equation for 


An 2 and Afl,j. 


(-R(l,3) - 
(-R(l,3) - K a | 2 ) 2 
• 


An 2 

• 

(-R«,3) - K A 2 ) k 




and AHj = An 2 


(n 3 + An 1 (R(2,3) + K A 
(n 3 + An 1 (R(2,3) + K A j)) 2 
• 

(n 3 + An 1 (R(2,3) + K A>3 )) k 

(Eq. D-106) 


g. If RMMAX(I) > YMMAX(I), solve following equations for A£j and 


(R(3,3) - 
(R(3,3) - K a ’ 1 ) 2 
• 

• 



- 

Os 

O 3 

+ A1 2 (R(2,3)+K JJj 
+ AH 2 (R(2,3) + K a | 3 T) 2 
• 

• 

• 

(R(3,3) - K A>1 ) k 




C + ^3 

• 

+ AJl 2 (R(2,3)+ K A ^ k 


(Eq. D-107) 


and An 2 = A l j 



(Eq. D-108) 
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4. If YMMAX(I) > PMMAX(I) and > RMMAX(I) 

a. Solve following matrix equation using a least square technique for 
Anj and An 2 


-(R(2,3) + K a 3 )j 
-<R(2,3) ♦ K a ’ 3 ) 2 
* 

(~R(1,3) - K a>2 )1 
(-R(1,3)- , K a , 2 )2 
• 


An^ 

— 

( n 3>l 

• 

• 

-<R<2,3)*K Ai3 ) k 

• 

• 

<-R(l,3) - K a 2 ) k 


i" 2 


• 

• 


An, 


b. x = + 


AY 


c. y = 


An 2 

AA 


(Eq. D-109) 
(Eq. D-110) 

(Eq. D-lll) 


d. A 1 1 = A n 2 


(Eq. D-112) 
(Eq. D-113) 

f. If PMMAX(I) > RMMAX(I), solve following equations for A m 9 

u 

and AH„ 


e. A = An^ 


(-R(i,3) - K A 2)l 




jjn 3 - Am 1 (K A}1 -R(3,3)7 

1 

<-R<1.3) - k a,2>2 
• 

• 


A m 2 

= 

(m 3 - AmjOC^j - R(3,3)T 
• 

• 

' 2 

• 

(-R(l,3) - K A 2 ) k 




• 

&3 - A mi (K A1 - R(3 ,3)j) 

k 


(Eq. D-114) 
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and AH j = A m 2 


g. If RMMAXG) > PMMAXO), solve the following equations for 


(R(2,3)- K a 3 ) t 
(R(2,3) - K a ’ 3 ) 2 


(R(2 ’ 3> - K A,3>k 


LI 2 anc * A m 2 





r — 





[+£, 3 - JL j(R(3,3) - K a j) 

1 




^H 3 -H 1 (R(3,3)-K A ’ i r 

2 


ah 2 

— 

• 

* 





• 

[a 3 -t l( R(3,3)-K Aii r 

|k 


(Eq. D- 115) 


and A m 2 = AH 2 



(Eq. D-116) 


5. If RMMAXO) > PMMAXfl) and > YMMAXfl) 


a. Solve the following matrix equation using a least squares technique 
for AH j and A 1 2 » 


(R(3,3)-K A1 ), 

«(3,3)-K a>1 ) 2 

• 

-(R (2,3) + K a Jj 
-(R( 2,3) + K a | 3 ) 2 
• 


Al, 


+(H 3 ) 1 

+(h 3 ) 2 

• 

• 

• 

<R<3.3> - K Ai ,) k 

• 

• 

-(R(2,3) + K A>3 ) k 


il 2 


• 

• 

+(H 3 ) k 


(Eq. D-117) 
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A A 1 


b. 

y= Air 

(Eq. D-118) 


ajl 2 


c. 

2 = AY 

(Eq. D-119) 

d. 

An 2 = AAj 

(Eq. D-120) 

e. 

A m 2 — A H 2 

(Eq. D-121) 


f. If PMMAX(I) > YMMAX(I), solve following equations for A 
and Anj. 


(K a>1 - R(3,3))j 
(K A,1 ' R(3,3)) 2 


Anij 


(m 3 - Am 2 (— R(l,3) - K A J) 

(m 3 - Am 2 (-R(l,3) - K a ’2 


1 

2 


(k a,i - R < 3 ’ 3 » k 


and An^ = Am^ 


(m 3 -4m 2 (-R(l,3)-K A J k 

(Eq. D-122) 


g. If YMMAX(I) > PMMAX(I), solve following equations for An^ and 

Anrij. 


-(R(2,3) t K Ai3>1 

- <R<2 ’ 3,+ K A,3 ) 1 
• 

• 


Anj 

= 

C"3 - An 2 '-“O' 3 ) - K A,2>) 1 

C"3 ' 4n 2 <- R(1 - 3 > - k aJ 2 

-(R(2, 3) + K A>3 ) k 




C"3 - 4n 2 <- R(1 ' 3) - K A,23 k 


(Eq. D-123) 
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and Am j = An^ 


h. x = 


Am j 
AN" 


(Eq. D-124) 


X. Calculation of Attitude Load Constants by Method II 

1. If FNMAXfl) > FYMAX(I), solve following equations for AN and AY. 


2 . 


*A,l ' R(3 - 3 »l 

<K A,I - r < 3 - 3 »2 
• 

• 


AN 

= 

(Mj); 

<N } ) ? 

• 

• 

• 

• 

(K a> , - R<3,3» k 




• 

• 


(Eq. D-125) 


and AY = AN 

If FYMAXO) > FNMAX, solve following equations for A Y and AN. 


-(R (2,3) + K A 3 )j 
-(R(2,3)+ K a ’ 3 ) 2 
• 

• 


AY 

= 

(Y 3>1 

(V 3>2 

• 

• 

• 

• 

-(R(2,3> . K A ,) k 




• 

• 

(Y 3>k 


(Eq. D-126) 


and A N = AY 
3. AA = AN-AW 


(Eq. D-127) 
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4. Determine Am^, A mg, An^, Ang, AJl^, AJlg, x, y, and z by calculation 
procedure given in Subsection W., item 3. 

Y. Balances Without Six Components 

1. For balances that do not have six components, set appropriate attitude 
tare constant to zero as indicated below. 

a. If balance does not have a normal-force component: AN = 0 

b. If balance does not have a axial-force component: A A = 0 

c. If balance does not have a side-force component: A Y = 0 

d. If balance does not have a pitching-moment component: 

Amj = Ariig = 0 

e. If balance does not have a rolling-moment component: 

AJlj = AJlg = 0 

f. If balance does not have a yawing moment component: 

An x = An 2 = 0 

Z. Initial Weight Tare Calculations 
1. Calculate initial weight tares 


a. 

N = 
0 

- 4N K A,1 

NFO 

(Eq. D-128) 

b. 

A = 
0 

- 4A K A,2 

AFO 

(Eq. D-129) 

c. 

m o ~ 

- Am l K A,l + 4nl 2 K A,2 

PMO 

(Eq. D-130) 


103 



D-52 


d. 

o 

II 

“l K A,l * At 2 K A,3 

RMO 

(Eq. D-131) 

e. 

n = 
0 

4n t K A,3 * 4n 2 K A,2 

YMO 

(Eq. D-132) 

f. 

o 

ii 

4YK A,2 

SFO 

(Eq. D-133) 


AA. New Values of Initial Weight Tares 

1. Go to Subsection V., item 5. and repeat calculation using new values of 
initial weight tares. Repeat iteration procedure until initial weight tares 
repeat to following accuracy. 


New - Old 
New 


< 0.005 


(Eq. D-134) 


BB. Point Calculations 

1. For each point, calculates 

a. = N 3 - [A N(K a>1 - R(3,3))J (Eq. D-135) 

b. = A 3 - [A A(-R(l,3) - K a>2 )J (Eq. D-136) 

c. = m 3 - j [Amj(K A)I - R(3,3))] - [a m 2 (R(l,3) + K^)] J 

’(Eq. D- 1 37) 

| 

d. VV |[4» 1 <R(3^-K Af ,»].[A» J (R(2 ( J).K Ai3 )]j 

* (Eq. D-138) 

'• W ![ 4 "l (R<2 - 3) * K A,3 ) ]-[ 4n 2 (R(1 ' 3) * K A,2 , ]| 

(Eq. D-139) 

f * Y (f = Y 3 + 0 Y(R(2,3) + K a>3 )] (Eq. D-140) 

i 

I 
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THETAS (0 S ) is measured in the 
tunnel or gravity X-Z plane. 
This is the current capability 
of the 16' TT model support 
system (\p s = 4> s = 0). e s is 
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(b) Tunnel support to undeflected 
balance axes. 




© 

V 0# * e k * 0# * *k ■ 180# 

© 

K» k a 0# * e k * 180 °» \ * °° 

© 

4> k ■ 0% 6 k * 180° + K°, «fr k » 180* 
or 

<J» k - 180°, e k " - K # . - 0 # 

or 

6 k * K *« \ * 180# * *k * 0# 

(c) Illustrations of knuckle angles. 
Figure D-l. Continued. 
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(d) Undeflected balance to balance axes. 
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(f) Illustration of primary balance to undeflected 
secondary balance rotations. 


Figure D-1. Continued. 
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(b) Case 1A, Case 1 held by opposite end. 




(c) Case 2, Balance rolled 90 (clockwise). 











(e) Case 3, Balance rolled 180 (inverted). 
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Figure D-2. Continued. 





( 9 ) Case 4, Balance rolled 90 (counterclockwise). 
Figure D-2. Continued. 
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(h) Case 4A, Case 4 held by opposite end. 
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(j) Case 5A, Case 5 held by opposite end. 
Figure D-2. Continued. 
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(1) Case 6A, Case 6 held by opposite end. 
Figure D- 2. Continued. 









(afterbody) model 



(n) Case 7A, Same as Case 7 held by opposite end. 
Figure D-2. Continued. 
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(p) Case 8A, Case 8 held by opposite end 
Figure D-2. Concluded. 




Fiaure 03. Definition of gravity and wind axes showing 
positive directions and rotation angles for 
wind to gravity transformations. 


Note: This figure assumes body axis 

origin is at balance pitch center. 


D-77 


fsl ^ 



129 


Fiqure D-4. Definition of balance and body axes showing positive 
directions and rotation angles for balance to model 
(body) transformations. 
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(b) Translation. 




Figure D-5 . Definition of angle of attack 
and angle of sideslip. 
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Rake Total Pressure/Static Pressure Assignments E-7 
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Figures 

E-l Definition of thrust angles E- 20 

E-2 Definition of Euler angles and directions E-21 
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MODULE E INTERNAL DRAG (OR EXIT-FLOW DISTRIBUTIONS) 


E-l 


SYMBOL 

AEXIT1 

AEXIT2 

ARAKE(I) 

CAI 

CAI1 

CAI2 

CDI 

CDIS 

CLI 

CNI 

CYI 

FTMDOT1 

FTMDOT2 

FTPR1 

FTPR2 

INDX(I,J) 


IRAKE 


KPR(I) 


NOMENCLATURE 

Exit areas for duct 1. sq. in. Not required for IRAKE = 2 or 
3. 

Exit areas for duct 2. sq. in. Not required for IRAKE = 2 or 
3. 

Exit area assigned to each rake total pressure PROBE(I), sq. 
in. Not required for IRAKE = 2 or 3. 

Total internal axial force coefficient. 

Internal axial force coefficient for duct 1. 

Internal axial force coefficient for duct 2. 

Total internal drag coefficient in the wind axis. 

Total internal drag coefficient in the stability axis. 

Total internal lift coefficient. 

Total internal normal force coefficient. 

Total internal side force coefficient. 

Mass flow rate at exit of duct 1, slugs/sec. 

Mass flow rate at exit of duct 2, slugs/sec. 

Ratio of nozzle exit total pressure to free stream static 
pressure for duct 1. 

Ratio of nozzle exit total pressure to free stream static 
pressure for duct 2. 

Table of values used to assign rake total pressures to 
specific static pressures, where I = static pressure probes 
assigned to J = table position. Not required for IRAKE = 2 
or 3. 

RAKE code. 

= o, set CAI=CDIS=CDI=0.0 and skip module 5. 

= 1, computes internal drag. 

= 2, measures exit flow distribution only. 

= 3, obtains internal drag from a given table. 

Needed to correct for bad rake static pressure probes. Set 
to 0.0 or 1.0, where I = static pressure probe. Not required 
for IRAKE = 2 or 3. 


preceding page blank not 
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E-2 

SYMBOL 

MEXIT1 

MEXIT2 

MODOT1 

MODOT2 

M/M01 

M/M02 

NPR1 

NPR2 

NPTR1 

NPTR2 

PD1/PTO 

PD2/PTO 

PRAKE(I) 

PR/PTO(I) 

PSIN1 

PSIN2 

PTD1/PTO 

PTD2/PTO 

PTRAKE(I) 

PTR/PTO(I) 


NOMENCLATURE 

Average exit mach number for duct 1. 

Average exit mach number for duct 2. 

Mass flow rate based on free-stream conditions for duct 1, 
slugs/sec. 

Mass flow rate based on free-stream conditions for duct 2, 
slugs/sec. 

Mass flow ratio for duct 1. 

Mass flow ratio for duct 2. 

Number of static pressure probes on the rake for duct 1. 
Maximum of 10. Not required for IRAKE = 3. 

Number of static pressure probes on the rake for duct 2. 
Maximum of 10-NPR1. Not required for IRAKE = 3. 

Number of total pressure probes on the rake for duct 1. 
Maximum of 50. Not required for IRAKE = 3. 

Number of total pressure probes on the rake for duct 2. 
Maximum of 50-NPTR1. Not required for IRAKE = 3. 

Ratio of the average duct static pressure to free-stream 
total pressure for duct 1. 

Ratio of the average duct static pressure to free-stream 
total pressure for duct 2. 

Rake static pressure, where I = probe number. 

Ratio of rake static pressure to free-stream total pressure, 
where I = probe number. 

Thrust axis yaw angle (degrees) for duct 1. Not required for 
IRAKE = 2 or 3. 

Thrust axis yaw angle (degrees) for duct 2. Not required for 
IRAKE = 2 or 3. 

Ratio of the average duct total pressure to free-stream 
total pressure for duct 1. 

Ratio of the average duct total pressure to free-stream 
total pressure for duct 2. 

Rake total pressure, where I = probe number. 

Ratio of rake total pressure to free-stream total pressure, 
where I = probe number. 
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SYMBOL 

NOMENCLATURE 

SCAP1 

Inlet capture area for duct 1. sq. in. Not required for 
IRAKE = 2 or 3. 

SCAP2 

Inlet capture area for duct 2 sq. in. Not required for 
IRAKE = 2 or 3. 

THETAN1 

Thrust axis Euler pitch angle (degrees), with respect to body 
axis for duct 1. Not required for IRAKE = 2 or 3. 

THETAN2 

Thrust axis Euler pitch angle (degrees), with respect to body 
axis for duct 2. Not required for IRAKE = 2 or 3. 
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APPENDIX E 


E-5 


Module E 

Internal Drag (or Exit Flow Distributions) 


Required Constants 

The constants for internal drag calculations are given in the nomenclatures. All 
constants are initialized to a value of 0.0. 


1. IRAKE - Rake code 

where IRAKE = 0, Set CAI = CDIS = CDI = 0.0 and skip this 

module. 

IRAKE = 1, compute internal drag 

IRAKE = 2, measure exit flow distribution only 

IRAKE = 3, obtain internal drag from a given table 

(Eq. E-l) 


(Eq. E-2) 

2. SCAP1, SCAP2 - inlet capture area, where SCAP1 is for duct 1 and 

SCAP2 is for duct 2. Not required for IRAKE = 2 or 

3. 


NPTR1 

£ ARAKE(I) = total exit area for duct 1 

1 = 1 

NPTR2 

£ ARAKE(I) = total exit area for duct 2 

I = NPTR1 + 1 


3 . AEXIT1, AEXIT2 - exit areas for ducts 1 and 2, respectively. Not 

required for IRAKE = 2 or 3. 

4. PSIN1, PSIN2 - Thrust axis yaw angle, with respect to body axis, for 

ducts 1 and 2, respectively. Positive direction is 
shown on Figure E-l. Not required for IRAKE = 2 or 
3. 
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5. THETAN1, THETAN2 - Thrust axis Euler pitch angles, with respect to 

body axis, for ducts 1 and 2, respectively, deg. 
Positive direction is shown on Figure E-l. Figure 
E-l also gives relations to obtain the Euler angle 
if not known directly. Not required for 
IRAKE = 2 or 3. 

.2 

6. AREF - Model reference area used for coefficients, in . 

If Module B or C is used, this constant is already 
specified. Not required for IRAKE = 2 or 3. 

Test for Module E Computations 

IF IRAKE = 0, skip module E. 

IF IRAKE = 3, do section T only. 

Rake Total Pressure 


1. Rake total pressures are called PTRAKE(I). Note that provisions are 
made to survey two exits at one time; however probes are numbered 
consecutively (max. of 50). For example, probes in the first exit may be 
numbered 1 through 30; probes in the second exit must start with number 
31. Where I = probe number. 

2. The ratio of rake total pressure to free-stream total pressure is called 
PTR/PTO(I), where PTO is from module A. 

3. The constants required from the project engineer are NPTR1, and NPTR2. 

Calculate PTR/PTO(I) for I = 1, NPTR1 + NPTR2 

PTRAKE(I) (e e _ 3) 

PTO 


PTR/PTO(I) = 


D. Rake Static Pressures 


E-7 


1. Rake static pressures are called PRAKE(I). Comments C.l. above apply 
except that the maximum number of probes is 10. 

2. The ratio of rake static pressure to free-stream total pressure is called 
PR/PTO(I). PTO is from module A. 

3. The constants required from the project engineer are NPR1, and NPR2. 

If NPR1 = 0, skip this part. 

Calculate PR/PTO(I) for I = 1, NPR1 + NPR2 

PR/PTO(I) = P p£P (I) - (Eq. E-4) 


E. Rake Total Pressure/Static Pressure Assignments 

1. If internal drag is to be computed, the project engineer must assign 
specific total pressure measurement to each static pressure measurement. 
This is done by supplying a table of I for all J, where I - total pressure 
measurements or probes which correspond to a specific J = static pressure 
measurement or probe. 

2. For example: 


J = 1 

I = 1, 2, 3, 4 

J = 2 

1= 5, 9, 11 

J = 3 
• 

1= 6, 7, 8, 10 
• 

• 

• 

J = NPR1 + NPR2 

• 

• 

I(Max) = NPTR1 + NPTR2 


3. The constants required from the project engineer are from the I, J table. 


141 



E-8 

If IRAKE = 2, skip this section. 

F. Duct Flow Static-to-Total-Pressure Ratio 

1. The ratio of duct flow static pressure to duct flow total pressure is called 
PR/PTR(J,I), where J and I are the combinations supplied in section E 
above. For the example shown in E., values of PR/PTR(J,I) are obtained 
for: PR/PTR1,1 
PR/PTR1,2 
PR/PTR1,3 
PR/PTR1,4 
PR/PTR2,5 
PR/PTR2,9 
PR/PTR2,11 
PR/PTR3,6 
PR/PTR3,7 
etc. 


If IRAKE = 2, skip this section. 

Do the following calculation for J = 1, NPR1 + NPR2 
Do the following calculation for I = those values assigned 


PR/PTR(J,I) = 


PR/PTCKJ) 

PTR/PTCKI) 


(Eq. E-5) 


G. Correct for Supersonic Duct Mach Numbers 

1. Local duct Mach number is called MD(I). Where I = total pressure probe 
number on which local Mach number is based. 

If IRAKE = 2, skip this section. 

Do the following calculation for J = 1, NPR1 + NPR2 
Do the following calculation for I = those values assigned 
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E-9 


If PR/PTR(3,I) < .5283, calculate MD(I) using the Newton Raphson method with 
an initial assumption of MDQ) = 1.0001 and correct the total pressure ratio for 
normal shock. 

/ 


H. Compute Subsonic Duct Mach Numbers 

1. This calculation is made for those I, J combinations for which 

PR/PTR (J, I)> .5283. 

If IRAKE = 2, skip this section. 

Do the following calculation for 3=1, NPR1 + NPR2 
Do the following calculation for I = those values assigned 
If PR/PTR(3,I)> .5283, calculate MDG) 

*>» ■ V 5 * [sfw] 2/7 - 5 E - 8> 

I. Compute Average Duct Pressure Ratios 

1. The ratio of the average duct total pressure to free-stream total pressure 
is called PTD1/PTO for duct 1 and PTD2/PTO for duct 2. 

2. The ratio of the average duct static pressure to free-stream total 
pressure is called PD1/PTO for duct 1 and PD2/PTO for duct 2. 
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3. The constants required from the project engineer are ARAKE(I), KPR(I), 
NPTR1, NPTR2, NPR1, and NPR2 


NPTR1 

l ARAKE(I) (j>TR/PTO(lj) 

PTDl/PTO = |^4 r! 

l ARAKE(I) 

1 = 1 


(Eq. E-9) 


NPTR1 

If l ARAKE(I) = 0.0, then PD1/PTO = 1.0 
1 = 1 

NPR1 

l KPR(I) [PR/PTOd^ 

PDl/PTO = 

l KPR(I) 

1 = 1 


(Eq. E-10) 


NPR1 

l KPR(I) = 0.0, then PDl/PTO = 1.0 
1 = 1 


If NPTR2 = 0.0, skip equations E-ll and E-12. 
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E-ll 


NPTR2 

l ARAKE(I) (PTR/PTO(I)) 

PTD2/PTO = ^ - f|J R1 + 1 (Eq. E-ll) 

l ARAKE(I) 

I = NPTR1 + 1 


NPTR2 

If l ARAKE(I) = 0.0, then PD2/PTO = 1.0 

I = NPTR1 + 1 


NPR2 

l KPR(I) (PR/PTO(I)) 

PD2/PT0 = (Eq. E-12) 

l KPR(I) 

1= NPR1 + 1 


NPR2 

If l KPR(I) = 0, then PD2/PTO = 1.0 

I = NPR1 + 1 

J. Mass-Flow Rates 

1. Mass-flow rate at the duct exit is called FTMDOT1 for duct 1 and 
FTMDOT2 for duct 2. 

2. Mass-flow rate based on free-stream conditions is called MODOT1 for 
duct 1 and MODOT2 for duct 2. 

3. TTO, MACH, and PO come from the tunnel parameters, module A. 
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The constants required from the project engineer are ARAKE(I)i SCAP1, 
SCAP2, NPTR1, and NPTR2. 

If IRAKE = 2, skip equations E-13, E- 14, E-15 and E-16. 


FTMDOT1 = 


NPTR1 

.028563 _ * l ARAKEU) * PRAKEO) * 

■^TTO + 459.67 I = 1 


PTR/PTO(I) 

PR/PTO0) 


1/7 


* MD(I) 


(Eq. E-13) 


where J corresponds to I from E.l. above. 


MODOT1 = 


(.028563) * (SCAP1) » (MACH) * (PO) * 
^TTO + 459.67 


1 + .2(MACH)" 


,1 I / 2 


(Eq. E-14) 


If NPTR2 = 0, skip equations E-15 and E-16. 


FTMDOT2 = 


NPTR2 

.028563 _ , l ARAKEQ) * PRAKEO) * 

^TTO + 459.67 I = NPTR1 + 1 


PTR/PTOfl) 

PR/PTOO) 


1/7 

* MD(I) 


(Eq. E-15) 


where J corresponds to I from E.l. above. 

11AfVvr , (.028563) » (SCAP2) ♦ (MACH) * (PO) * 
MODOT2 = — : 

^TTO + 459.6 7 


1 + .2(MACH) 2 


1/2 


(Eq. E-16) 


K. Mass-Flow Ratio 

1. Mass-flow ratio for duct 1 is called M/M01. Mass-flow ratio for duct 2 is 
called M/M02.. 

If IRAKE = 2, skip the remainder of this section. 
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M/MOi = 


FTMDOT1 

MODOT1 


If MODOT1 = 0.0, M/M01 = 0.0 


If NPTR2 = 0, skip equation E-18. 


M/MO o FTMDOT2 
M/M02 = modot2 


L. Free-Stream Velocity 

1. Free-stream velocity is called VO. 

2. TTO and MACH are from module A. 

VQ _ 49.01428>/tTO + 456.67 

Vi + .2 (MACH) 2 


(Eq. E-17) 


(Eq. E-18) 


(Eq. E-19) 


M. Average Exit Mach Number 

1. The average exit Mach number for duct 1 is always called MEXIT1. The 

average exit Mach number for duct 2 is always called MEXIT2. 


If IRAKE = 2, skip the remainder of this section. 


MEXIT1 = 



PTD1/PTO 

pdi/ptcT 


2/7 


-5 


(Eq. E-20) 


If NPTR2 = 0, skip equation E-21. 


MEXIT1 = 


J c 1 

PTD2/PTO 

2/7 

P * 

PD2/PTO 

— J 


(Eq. E-21) 


147 


Internal Axial Force 

1. The internal axial force is called All and AI2 for ducts 1 and 2, 
respectively. 

2. The internal axial force coefficient is called CAI1 and CAI2 for ducts 1 
and 2, respectively. 

3. The total internal axial force coefficient is called CAI. 

4. PTO, PO, and QO are from the tunnel parameters, module A. 

5. PSI and THETA are from the balance and weight tare calculations, module 
D. Positive directions for PSI and THETA are shown on Figure E-2. 

6. The constants required from the project engineer are AEXIT1, AEXIT2, 
PSINl, PSIN2, THETAN1, THETAN2, AREF, and NPTR2 


If IRAKE = 2, skip the remainder of this section. 


All = [(FTMDOT1) * VO * COS(PSI) * COS(THETA)] 

-|[l.4*(PDl/PTO) * PTO * (MEXIT1) 2 ] + [((PD1/PTO) * PTO) - Pojj 
* (AEXIT1) * COS(PSINl) * COSfTHETANl) (Eq. E-22) 


CAI1 = 


All 

(QO)*(AREF) 


(Eq. E-23) 


CAI = CAI 1 (Eq. E-24) 

If NPTR2 = 0, skip equations E-25, E-26 and E-27. 

AI2 = [(FTMDOT2) * VO * COS(PSI) * COStfHETA)] 

-|[l.4*(PD2/PTO) * PTO * (MEXIT2) 2 ] + [((PD2/PTO) * PTO) - Po]j 
* (AEXIT2) * COS(PSIN2) * COS(THETAN2) (Eq. E-25) 


E- 15 


CAI2 * (QO) « (AREF) C * E ' 26) 

CAI = CAI1 + CAI2 (Eq. E-27) 

O. Internal Normal Force 


1. The internal normal force is called Nil and NI2 for ducts 1 and 2, 
respectively. 

2. The internal normal force coefficient is called CNI1 and CNI2 for ducts 1 
and 2, respectively. 

3. The total internal normal force coefficient is called CNI. 

4. PTO, PO, and QO are from the tunnel parameters, module A. 

5. PSI, THETA, and PHI are from the balance and weight tare calculations, 
module D. Positive directions are shown on Figure E-2. 

6. The constants required from the project engineer are AEXIT1, AEXIT2, 
THETAN 1, THETAN2, AREF, and NPTR2 

If IRAKE = 2, skip the remainder of this section. 


Nil = 

| (FT M DOT 1)* VO* jjcoS(PHI)*SIN(THETA)*COS(P5I)+SIN(PHI)*SIN(PSI)] j 
+ |[l.4*(PDl/PTO)*PTO*(MEXITl) 2 ]<- [((PDl/PTO)*PTO) - po)| 

* (AEXIT 1)*SIN(THETAN1) (Eq. E-28) 

CN11 = (QO) ** (AREF) ® q * E_29) 


CNI = CNI1 (Eq. E-30) 

If NPTR2 = 0, skip equations E-31, E-32 and E-33. 
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NI2 = | (FT M DOT 2)* VO* [cOS(PHI)*SIN(THETA)*COS(PSI) + SIN(PHI)*SIN(PSlf] j 

+ j jj.4*(PD2/PTO)*PTO*(MEXIT2) 2 ]+ [((PD2/PTO)*PTO) - Po] j 
* (AEXIT2)*SIN(THETAN2) E ’ 31 ^ 


CNI2 = 


NI2 

(QO) * (AREF) 


(Eq. E-32) 


CNI = CNI1 + CNI2 


(Eq. E-33) 


P. Internal Side Force 

1. The internal side force is called YI1 and YI2 for ducts 1 and 2, 
respectively. 

2. The internal side force coefficient is called CYI1 and CYI2 for ducts 1 
and 2, respectively. 

3. The total internal side force coefficient is called CYI. 

it. PTO, PO, and QO are from the tunnel parameters, module A. 

5. PSI, THETA, and PHI are from the balance and weight tares calculations, 
module D. 

6. The constants required from the project engineer are AEXIT1, AEXIT2, 
THETAN 1, THETAN2, PSIN1, PSIN2, AREF, and NPTR2. 


If IRAKE = 2, skip the remainder of this section. 

Yll = | (FT M DOT 1)* VO* [siN(PHI)*SIN(THETA)*COS(PSI)-COS(PHI)*SIN(PSI)] ] 

+ j |1.4*(PDl/PTO)*PTO*(MEXITl) 2 ]+ [((PD 1 /PTO)* PTO) - Po] ] 

* (AEXIT l)*COS(THETANl)*SIN(PSIN 1) (Eq. E.34) 

r Yll - YI 1 (Eq. E-35) 

(QO) * (AREF) 
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CYI = CYI1 


E-17 

(Eq. E-36) 


If NPTR2 = 0, skip equations E-37, E-38 and E-39. 

YI2 = j (FTMDT2)*VO* |siN(PHI)*SIN(THETA)*COS(PSI)-COS(PHI)*SIN(PSI)] j 
+ | [i.4*(PD2/PTO)*PTO*MEXIT2) 2 ]+[((PD2/PTO)*PTO) - Po] j 
* (AE XIT 2 )*COS(THETAN 2 )*SIN (PSIN 2) (Eq* E-37) 


CYI2 = 


YI2 

(QO) * (A REF) 


(Eq. E-38) 


CYI = CYI1 + CYI2 


(Eq. E-39) 


Q. Flow-Through Pressure Ratio 

1. The nozzle exit (flow-through) total pressure in ratio to free-stream 
static pressure is called PTD1/PO and PTD2/PO for ducts 1 and 2, 
respectively. 

2. PTO and PO are from tunnel parameters, module A. 

If IRAKE = 2, skip the remainder of this section. 

PTDl/PO . (PTDI/ p »JO).(PTO L Eq. E-*0> 

If NPTR2 = 0, skip equation E-41. 

PTD2/PO = (PTD2/PTO) * (FTP) , <Eq.E-U) 


R. Internal Drag 

1. The internal drag coefficient based on the stability axes is called CDIS1 
and CDIS2 for ducts 1 and 2, respectively. 
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E-18 


2 . 


The total internal drag coefficient in the stability axis is called CDIS. 


3. ALPHA is from the balance and weight tares computations, module D. 

4. The internal drag coefficient based on the wind axis is called CDI1 and 
CDI2 for ducts 1 and 2, respectively. 

5. The total internal drag coefficient in the wind axis is called CDI. 

6. BETA is from module D. 

If IRAKE = 2, skip the remainder of this section. 

CDIS1 = (CNI1) * SIN(ALPHA) + (CAI1) * COS(ALPHA) (Eq. E-42) 

CDI1 = (CDIS1) * COS(BETA) - (CYI1) * SIN(BETA) (Eq. E-43) 

CDIS = CDIS1 (Eq. E-44) 

CDI = CDU (Eq. E-45) 

If NPTR2 = 0, skip equations E-46, E-47, E-48 and E-49. 

CDIS2 = (CNI2) * SIN(ALPHA) + (CAI2) * COS(ALPHA) (Eq. E-46) 

CDI2 = (CDIS2) * COS(BETA) - (CYI2) * SIN(BETA) (Eq. E-47) 

CDIS = CDIS1 + CDIS 2 (Eq. E-48) 

CDI = CDU + CDI2 (Eq. E-49) 
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Internal Lift 


E-19 


S. 


1. The internal lift coefficient based on stability axis (also wind axis) is 
called CLI1 and CLI2 for ducts 1 and 2, respectively. 

2. The total internal lift coefficient is called CLI. 

3. ALPHA is from the balance and weight tares computations in module D. 

If IRAKE = 2, skip the remainder of this section. 

CLI1 = (CNI1) * COS(ALPHA) - (CAI1) * SIN(ALPHA) (Eq. E-50) 

CLI = CLI1 (Eq. E-51) 


If NPTR2 = 0, skip equations E-52 and E-53. 


CLI2 = (CNI2) * COS(ALPHA) - (CAI2) * SIN(ALPHA) (Eq. E-52) 

CLI = CLI1 + CLI2 (Eq. E-53) 


T. Internal Drag and Axial Force Tables 
If IRAKE ^ 3, skip this section. 

CAI, CDSI, and CDI are supplied in tables as functions of MACH, ALPHA and 
PSI. 
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plane of paper, but 
rotated about line BO] 



Figure E-2. Definition of Euler angles and directions. 
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APPENDIX F 

Pressure Coefficients and Integrated Forces 


Nomenclatures F-l 

Required Constants F-5 

Test For Module F Computations F-5 

Free-Stream Static and Dynamic Pressures F-5 

Coefficient Calculations F-6 

Total Pressure Drag Coefficient F-6 

Internal Static Pressure Ratio F-7 
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MODULE F PRESSURE COEFFICIENT AND INTEGRATED FORCES 


F-l 


SYMBOL 

ARAAU(I) 

ARABU(I) 

ARADU(I) 

ARAFU(I) 

ARAGU(I) 

ARAHU(I) 

ARANU(I) 

ARASU(I) 

AREF 

ARHAU(I) 

ARHBU(I) 

ARHDU(I) 

ARHFU(I) 

ARHGU(I) 

ARHHU(I) 

ARHNU(I) 

ARHSU(I) 

ARNAU(I) 

ARNBU(I) 

ARNDU(I) 

ARNFU(I) 

ARNGU(I) 

ARNHU(I) 

ARNNU(I) 

ARNSU(I) 

ARPAU(I) 

ARPBU(I) 

ARPDU(I) 

ARPFU(I) 

ARPGU(I) 

ARPHU(I) 

ARPNU(I) 

ARPSU(I) 


Axial 

Force 


Hinge 

Moment 


Normal 

Force 


Pitch 

Moment 


NOMENCLATURE 


Areas to be used with pressure groups to compute 
integrated forces, where I = orifice number, sq. in. 


Model reference area from module B, sq. in. 


Areas to be used with pressure groups to compute 
integrated forces, where I = orifice number, sq. in. 


Areas to be used with pressure groups to compute 
integrated forces, where I = orifice number, sq. in. 


Area times moment arm, sq. in. to be used with pressure 
group to compute integrated moments, I = orifice, sq. in. 


PSfcCEDWQ FA3E SUNMk NOT 
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F-2 

SYMBOL 

CBAR 

CDAUN 

CDBUN 

CDDUN 

CDFUN 

CDGUN 

CDHUN 

CDNUN 

CDSUN 

CDPR 

CFAUN 

CFBUN 

CFDUN 

CFFUN 

CFGUN 

CFHUN 

CFNUN 

CFSUN 

CHMAUN 

CHMBUN 

CHMDUN 

CHMFUN 

CHMGUN 

CHMHUN 

CHMNUN 

CHMSUN 

CLAUN 

CLBUN 

CLDUN 

CLFUN 

CLGUN 

CLHUN 

CLNUN 


NOMENCLATURE 

Pitching moment reference length from module D, in. 
Integrated pressure drag coefficients. 


Total integrated drag coefficient. 

Integrated pressure axial force coefficients. 


Integrated pressure hinge moment coefficients. 


Integrated pressure lift coefficients. 


▼ 
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F-3 


SYMBOL 
CLSUN 
CL PR 
CNAUN 
CNBUN 
CNDUN 
CNFUN 
CNGUN 
CNHUN 
CNNUN 
CNSUN 
CPMAUN 
CPMBUN 
CPMDUN 
CPMFUN 
CPMGUN 
CPMHUN 
CPMNUN 
CPMSUN 
CPMPR 
KCDA 
KCDB 
KCDD 
KCDF 
KCDG 
KCDH 
KCDN 
KCDS 
PAUN 
PBUN 
PDUN 
PFUN 
PGUN 
PHUN 
PNUN 
PSUN 



NOMENCLATURE 


Total integrated lift coefficient. 

Integrated pressure normal force coefficients. 


Integrated pressure pitching moment coefficients. 


Total integrated pitching moment coefficient. 


Constants provided by the engineer. 

(0.0 

or 

1.0) 

Constants provided by the engineer. 

(0.0 

or 

1.0) 

Constants provided by the engineer. 

(0.0 

or 

1.0) 

Constants provided by the engineer. 

(0.0 

or 

1.0) 

Constants provided by the engineer. 

(0.0 

or 

1.0) 

Constants provided by the engineer. 

(0.0 

or 

1.0) 

Constants provided by the engineer. 

(0.0 

or 

1.0) 

Constants provided by the engineer. 

(0.0 

or 

1.0) 


Individual pressures to be 

used with each type of pressure coefficient 

for computation of integrated forces and 

moments. Maximum number of each type is 125. lbs/sq. in. 
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F-4 

SYMBOL NOMENCLATURE 

PRATI(II) Ratio of nozzle internal static pressure to nozzle total 

pressure, where II = orifice number (PGUN only). 
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APPENDIX F 


F-5 


Module F 

Pressure Coefficients and Integrated Forces 

Eight groups of pressure coefficients may be computed under this module. Names assigned 
to each group are arbitrary. Final names may be inserted with finalized data printout 
headers. These groups are PAUN, PBUN, PDUN, PFUN, PGUN, PHUN, PNUN, and PSUN. 

A. Required Constants 

The required constants for module F are given in the nomenclatures. 

1. All constants are initialized to a value of zero. The project engineer need 
only supply those constants which are required for those quantities to be 
computed. 

2. KAUN, KBUN, KDUN, KFUN, KGUN, KHUN, KNUN, KSUN - number of 
individual pressures to be used with each type of pressure coefficient for 
computation of integrated forces and moments. (125 maximum for each 
type) 

B. Test for Module F Computations 
If KPRESS = 0, skip this module. 

C. Free-Stream Static and Dynamic Pressures 

Free-stream static and dynamic pressures to be used for computing pressure 
coefficients are obtained from module A; however, for individual pressure 
transducers, an average value is used. For Scanivalve pressure transducers, 
PO and QO are calculated for each frame of data (one frame/port). 
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D. Coefficient Calculations 


CPSUN(I) = (PSUN(I) - PO)/QO (Eq. F-l) 

KSUN 

CFSUN = l CPSUN(I) * ARASU(I)/AREF (Eq. F-2) 

1= 1 

KSUN 

CNSUN = l CPSUN(I) * ARNSU(I)/AREF (Eq. F-3) 

I = 1 

KSUN 

CP MS UN = l CPSUN(I) * ARPSU(I)/AREF * CBAR (Eq. F-4) 

I = 1 

KSUN 

CHMSUN = l CPSUN(I) * ARHSU(I)/AREF * CBAR (Eq. F-5) 

I = 1 

CDSUN = CFSUN * COS (ALPHA) + CNSUN * SIN (ALPHA) (Eq. F-6) 

CLSUN = CNSUN * COS (ALPHA) - CFSUN * SIN (ALPHA) (Eq. F-7) 

These equations are the same for all pressure groups. 


E. Total Pressure Drag Coefficient 

CDPRessure = KCDS * CDSUN + KCDA * CDAUN + KCDB * CDBUN + KCDN * 
CDNUN + KCDD * CDDUN + KCDF * CDFUN + KCDG * CDGUN + KCDH * 
CDHUN 

(Eq. F-8) 

where KCDS, KCDB, KCDA, KCDN, KCDD, KCDF, KCDG, KCDH are constant 
inputs, either 0 or 1.0. 


F-7 

F. Internal Static Pressure Ratio 

PRATI(II) = PGUN(II)/PTENG1 (Eq. F-9) 

NOTE : In addition to the pressure coefficients, this ratio is for PGUN 

measurements only. 
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APPENDIX G 
Thrust Removal Options 


Nomenclatures G-l 

General Information G-5 

Required Constants G-5 

Quantities Required G-5 

Compute Thrust and Static Thrust Terms 
IF = 1 G-6 

Single Balance/All Metric IF1 = 1 G-ll 

Single Balance/Afterbody Metric IF2 = 1 G-l 3 

Two-Balance/Afterbody Metric IFAF1 = 1 G-15 

Two-Balance/Afterbody Metric IFAF2 = 1 G-18 

Two-Balance/Afterbody Metric IFAFN1 = 1 G-20 

Two-Balance/Afterbody Metric IFAFN2 = 1 G-22 

Single Balance, Thrust Removal All 
Components IFAF G-24 

When IFAFN = 1 G-28 

Bifurcate Support Mode Two Balance/Afterbody 
Metric IDN = 1 G-29 

Other Options G-33 


8UWK NOT 
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MODULE G THRUST REMOVAL OPTIONS 


G-l 


SYMBOL 

AEX 

CAAERO 

CASCADE 

CDAERO 

CDNOZ 

C(F-ANOZ) 

C(F-DNOZ) 

CDSAER 

CDWAER 

CF 

CF/CFI 

CFJ 

CFJC 

CFJET 

C(F-A) 

C(F-D) 

CLAERO 

CLJET 

CLNOZ 

CLNOZT 

CLSAER 

CLWAER 

CMAERO 

CMJ 

CMJC 

CM JET 

CMNOZ 

CMNOZT 

CMSAER 

CMWAER 

CNAERO 

CNJ 

CNJC 


NOMENCLATURE 
Nozzle exit area. 

Thrust removed axial force coefficient. 

Resultant angle of jet exhaust, degrees. 

Thrust removed drag coefficient. 

Nozzle drag. 

Thrust minus nozzle axial force coefficient. 

Thrust minus nozzle drag coefficient. 

Thrust removed stability axis drag coefficient. 

Thrust removed wind axis drag coefficient. 

Jet axial force coefficient (from balance and pressures). 
Ratio of thrust (from balance and pressures) to ideal thrust. 
Jet axial force coefficient. 

Computed jet axial force coefficient. 

Jet reaction axial force coefficient. 

Thrust minus axial force coefficient. 

Thrust minus drag coefficient. 

Thrust removed lift coefficient. 

Jet reaction lift coefficient. 

Thrust removed nozzle lift. 

Nozzle lift plus thrust. 

Thrust removed stability axis lift coefficient. 

Thrust removed wind axis lift coefficient. 

Thrust removed pitching moment coefficient. 

Jet pitching moment coefficient. 

Computed jet pitching moment coefficient. 

Jet reaction pitching moment coefficient. 

Thrust removed nozzle pitching moment. 

Nozzle pitching moment plus lift. 

Thrust removed stability axis pitching moment coefficient. 
Thrust removed wind axis pitching moment coefficient. 
Thrust removed normal force coefficient. 

Jet normal force coefficient. 

Computed jet normal force coefficient. 


PR* CEE»«G BLANK NOT FUMED 
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G-2 

SYMBOL 

CRAERO 

CRJC 

CRJET 

CRSAER 

CRWAER 

CSAERO 

CSJC 

CSJET 

CSSAER 

CSWAER 

CT 

CTS 

CTST 

CTSY 

CTS2 

CTS2T 

CTS2Y 

CTT 

CTY 

CYAERO 

CYJC 

CYJET 

CYSAER 

CYWAER 

DELTA 

DELTAY 

DELTA1 

DELTA2 


NOMENCLATURE 

Thrust removed rolling moment coefficient. 

Computed jet rolling moment coefficient. 

Jet reaction rolling moment coefficient. 

Thrust removed stability axis rolling moment coefficient. 
Thrust removed wind axis rolling moment coefficient. 

Thrust removed side force coefficient. 

Computed jet side force coefficient. 

Jet reaction side force coefficient. 

Thrust removed stability axis side force coefficient. 

Thrust removed wind axis side force coefficient. 

Computed resultant thrust coefficient about pitch axis. 
Resultant static thrust coefficient, main balance, about 
pitch axis. 

Resultant static thrust coefficient, main balance. 

Resultant static thrust coefficient, main balance, about yaw 
axis. 

Resultant static thrust coefficient, second balance, about 
pitch axis. 

Resultant static thrust coefficient, second balance. 

Resultant static thrust coefficient, second balance, about 
yaw axis. 

Computed resultant thrust coefficient. 

Computed resultant thrust coefficient about yaw axis. 

Thrust removed yawing moment coefficient. 

Computed jet yawing moment coefficient. 

Jet reaction yawing moment coefficient. 

Thrust removed stability axis yawing moment coefficient. 
Thrust removed wind axis yawing moment coefficient. 
Computed thrust vector angle about pitch axis, degrees. 
Computed thrust vector angle about yaw axis, degrees. 

Static thrust vector angle, main balance, about pitch axis, 
degrees. 

Static thrust vector angle, second balance, about pitch axis, 
degrees. 


170 


; « 




G-3 


SYMBOL 

DELT1Y 

DELT2Y 

ETAABS 

(F-A)/FI 

(F-ANOZ)/FI 

FGT/FI 

FGT2/FI 

FGY/FI 

FG/FI 

FG2/FI 

FG2Y/FI 

FJ1/FI 

FJ2/FI 

(F-D)/FI 

(F-DNOZ)/FI 

FN/FI 

FT/FI 

F/FI 

IDA 

IDN 

IF 

IF1 

IF2 

IFAF 

IFAF1 

IFAF2 

IFAFN 

IFAFN1 

IFAFN2 

LENGTH(I) 


NOMENCLATURE 

Static thrust vector angle, main balance, about yaw axis, 
degrees. 

Static thrust vector angle, second balance, about yaw axis, 
degrees. 

Isentropic vacuum or stream thrust coefficient. 

Ratio of thrust minus axial force to ideal thrust. 

Ratio of thrust minus nozzle axial force to ideal thrust. 

Total static resultant thrust ratio, main balance. 

Total static resultant thrust ratio, second balance. 

Static resultant thrust ratio, main balance, about yaw axis. 
Static resultant thrust ratio, main balance, about pitch axis. 
Static resultant thrust ratio, second balance, about pitch 
axis. 

Static resultant thrust ratio, second balance, about yaw 
axis. 

Static thrust ratio, main balance. 

Static thrust ratio, second balance. 

Ratio of thrust minus drag to ideal thrust. 

Ratio of thrust minus nozzle drag to ideal thrust. 

Ratio of normal force to ideal thrust. 

Total resultant thrust ratio. 

Ratio of thrust to ideal thrust. 

Engineer's option. 

Future option. 

Computes thrust and static thrust terms when IF=1. 
Computes single balance/all metric when IF1=1. 

Computes single balance/afterbody metric when IF2=1. 

Single balance, thrust removal from all components. 
Computes two balances/afterbody metric when IFAF1=1. 
Computes two balances/afterbody metric when IFAF2=1. 
Future option. 

Computes two balances/afterbody metric when IFAFN1=1. 
Computes two balances/afterbody metric when IFAFN2=1. 
Lengths to transfer moments to relative station. 
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G-4 

SYMBOL 

PM/FI 

RM/FI 

SF/FI 

SPLAY 

SPLAY1 

YM/FI 


NOMENCLATURE 

Ratio of pitching moment to ideal thrust. 
Ratio of rolling moment to ideal thrust. 
Ratio of side force to ideal thrust. 
Projected roll angle of jet exhaust, degrees. 
Projected roll angle of jet exhaust, degrees. 
Ratio of yawing moment to ideal thrust. 
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G-5 

APPENDIX G 
Module G 

Thrust Removal Options 


A. General Information 

The following options are used to remove thrust and to obtain various 
aerodynamic and aeropropulsion parameters usually required for most 16-Ft. 
Transonic Tunnel investigations. The various constants are keyed to typical 
balance arrangements used and may be used for most test setups. This section 
requires computed inputs from modules A, B, C, D and E. The engineer should 
refer to each module for exact definition of the computed quantity. These 
options will work for both fully and partially metric models for both 
longitudinal and lateral data. 

B. Required Constants 

1. IF, IF 1 , IF2, IFAF, IFAF1, IFAF2, IFAFN, IFAFN1, IFAFN2, ID and IDN. 

C. Quantities Required 

1. MODULE A 

a. PO & QO 

2. MODULE B 

a. NPR 

b. CFI 

3. MODULE C 

a. CDF AFT - afterbody + nozzle skin friction 

b. CDFNOZ - nozzle skin friction 


173 


4 


MODULE D 


1. ALPHA 

2. CN1, CA1, CMY1 \ 

CY1, CMX1, CMZ1 > MAIN BALANCE 

3. CDS1, CLS1 ' 

4. CN2, CA2, CMY2 ^ 

CY2, CMX2, CMZ2 > SECOND BALANCE (2) 

5. CDS2, CLS2 ’ 


5. MODULE F 

1. CFSUN, CNSUN, CPMSUN 

2. CDSUN, CLSUN 

3. CFBUN, CNBUN, CPMBUN 

4. CDBUN, CLBUN 


} 

} 


AFTERBODY PRESSURE FORCES 


NOZZLE PRESSURE FORCES 


Compute Thrust and Static Thrust Terms IF = 1 


1. Compute Thrust 

a. If NPR < 1.2, CFJ = CNJ = CMJ = 0 = CRJ = CYJ = CSJ 

b. The computed jet axial force coefficient is 

CFJC = ^ * (kCFJ(NPR) + ICFj) (Eq. G-l) 

c. The computed jet normal force coefficient is 

CNJC = * (KCNJ(NPR) + ICNJ] (Eq. G-2) 


d. 


The computed jet pitching moment coefficient is 



G-7 

CMJC = ^ * (KCMJ(NPR) + ICMj) (Eq. G-3) 

e. The computed jet rolling moment coefficient is 

CRJC = ^ * (KCRJ(NPR) + ICRj) (Eq. G-4) 

f. The computed jet yawing moment coefficient is 

CYJC = * (KCYJ(NPR) + ICYJ3 (Eq. G-5) 

g. The computed jet side force coefficient is 

CSJC = * (kcSJ(NPR) + ICSJ3 (Eq. G-6) 

h. Table input is as follows: (Need six tables) 

Up to five values per table may be used. 

NPR Range Slope Intercept 

Compute Static Thrust Terms 

a. The resultant static thrust coefficient about the pitch axis for the 
main balance is 

CTS = Jem 2 + CA1 2 (Eq. G-7) 

b. The resultant static thrust ratio about the pitch axis for the main 
balance is 

FG/FI = CTS/CFI (Eq. G-8) 
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G-8 


c. The static thrust ratio for the main balance is 

FJ1/FI = -CA1/CFI (Eq. G-9) 

d. The static thrust vector angle about the pitch axis for the main 
balance is 

DELTA1 = TAN" 1 (-CN1/CA1) (Eq. G-10) 

e. The resultant static thrust coefficient about the yaw axis for the 
main balance is 

CTSY = v/cYl 2 + CA1 2 (Eq. G-ll) 

f. The resultant static thrust ratio about the yaw axis for the main 
balance is 

F GY/FI = CTSY/CFI (Eq. G-12) 

g. The static thrust vector angle about the yaw axis for the main 
balance is 

DELTA1Y = TAN -1 (-CY1/CA1) (Eq. G-13) 

h. The resultant static thrust coefficient for the main balance is 

CTST = / CN1 2 + CA1 2 + CY1 2 (Eq. G-14) 

i. The total resultant static thrust ratio for the main balance is 

FGT/FI = CTST/CFI (Eq. G-15) 
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j- 


The isentropic vacuum thrust or stream thrust coefficient is 


computed by 


ETAABS = 


FGT/FI * FI AEX 

PTJAVG PTJ/PO 

PE/PTJ * AEX * (1 + Y * ME) 2 


(Eq. G-16) 


where AS = WPWITO * AT(1) 


k. The nozzle exit Mach number ME is computed from 


AS/AEX 


= HI ME^l + 0.2 ME 2 ^ 


-3 


and the nozzle exit pressure ratio is from 

( 2V 7 /2 

PE/PTJ =\1 + 0.2 ME / 


(Eq. G-17) 


1. The resultant static thrust coefficient about the pitch axis for the 
second balance is 

CTS2 = v/cN2 2 + CA2 2 (Eq. G-18) 


m. The resultant static thrust ratio about the pitch axis for the second 
balance is 


FG2/FI = CTS2/CFI (Eq. G-19) 

n. The static thrust ratio for the second balance is 

FJ2/FI = -CA2/CFI (Eq. G-20) 

o. The static thrust vector angle about the pitch axis for the second 
balance is 

DELTA2 = TAN -1 (-CN2/CA2) (Eq. G-21) 
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p. The resultant static thrust coefficient about the yaw axis for the 
second balance is 


CTS2Y = 


/ 


CY2 2 + 


CA2' 


(Eq. G-22) 


q. The resultant static thrust ratio about the yaw axis for the second 
balance is 


FG2Y/FI = CTS2Y/CFI (Eq. G-23) 

r. The static thrust vector angle about the yaw axis for the second 
balance is 

DELTA2Y = TAN -1 (-CY2/CA2) (Eq. G-24) 

s. The resultant static thrust coefficient for the second balance is 

CTS2T = J CN2 2 + CA2 2 + CY2 2 (Eq. G-25) 

t. The total resultant static thrust ratio for the second balance is 


FGT2 = CTS2T/CFI 


(Eq. G-26) 


u. The splay angle is 


SPLAY = ATAN(CY1/CN1) 


(Eq. G-27) 


v. The cascade angle is 


CASCADE = TAN 


TAN (DELTA1) 


\J\ + TAN 2 (SPLAY) * TAN 2 (DELTA1) 


(Eq. G-28) 
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G-ll 


w. The ratio of normal force to ideal thrust is 

FN/FI = CN1/CFI (Eq. G-29) 

x. The ratio of side force to ideal thrust is 

SF/FI = CY1/CFI (Eq. G-30) 

y. The ratio of rolling moment to ideal thrust is 

RM/FI = CMX1/(CFI * LENGTH1) (Eq. G-31) 

z. The ratio of pitching moment to ideal thrust is 

PM/FI = CMY1/(CFI * LENGTH2) (Eq. G-32) 

aa. The ratio of yawing moment to ideal thrust is 
YM/FI = CMZ1/(CFI * LENGTH3) 

(Eq. G-33) 


E. Single Balance/ All Metric IF1 = 1 
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The resultant thrust coefficient about the pitch axis is 


CT = \/cNJC 2 + CFJC 2 

The thrust vector about the pitch axis is 

DELTA = TAN' 1 (CNJC/CFJC) 

The jet reaction lift coefficient is 
CLJET = CT(SIN (ALPHA + DELTA]) 

The jet reaction axial force coefficient is 
CFJET = CT (COS (ALPHA + DELTA]) 

The thrust removed lift coefficient is 
CLAERO = CLS1 - CLJET 
The thrust removed drag coefficient is 
CDAERO = CDS1 + CFJET 

The thrust removed pitching moment coefficient is 

CMAERO = CMYS1 - CMJC 

The thrust minus axial force coefficient is 

C(F - A) = -CA1 

The thrust minus drag coefficient is 


(Eq. G-34) 


(Eq. G-35) 


(Eq. G-36) 


(Eq. G-37) 


(Eq. G-38) 


(Eq. G-39) 


(Eq. G-40) 


(Eq. G-41) 


C(F - D) = -CDS1 


(Eq. G-42) 


10. The ratio of thrust minus axial force to ideal thrust is 
(F - A)/FI = C(F - A)/CFI 

11. The ratio of thrust minus drag to ideal thrust is 
(F - D)/FI = C(F - D)/CFI 

12. The ratio of thrust to ideal thrust is 
F/FI = CFJC/CFI 

13. The thrust minus nozzle drag coefficient is 
C(F - DNOZ) = CFJET - CDBUN 

14. The ratio of thrust minus nozzle drag to ideal thrust is 
(F - DNOZ)/FI = C(F - DNOZ)/CFI 

F. Single Balance/Afterbody Metric IF2 = 1 


G-13 

(Eq. G-43) 
(Eq. G-44) 
(Eq. G-45) 


(Eq. G-46) 


(Eq. G-47) 
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G-14 


1. The resultant thrust coefficient about the pitch axis is 

CT = n/cNJC 2 + CFJC 2 (Same as Eq. G-34) 

2. The thrust vector angle about the pitch axis is 

DELTA = TAN -1 (CNJC/CFJC) (Same as Eq. G-35) 

3. The jet reaction lift coefficient is 

CLJET = CTfsiN (ALPHA + DELTA^| (Same as Eq. G-36) 

4. The jet reaction axial force coefficient is 

CFJET = CT (COS (ALPHA + DELTA)) (Same as Eq. G-37) 

5. The thrust removed lift coefficient is 

CLAERO = CLS1 - CLJET (Same as Eq. G-38) 

6. The thrust removed drag coefficient is 

CDAERO = CDS1 + CFJET (Same as Eq. G-39) 

7. The thrust removed pitch moment coefficient is 

CMAERO = CMYS1 - CMJC (Same as Eq. G-40) 

8. The thrust minus axial force coefficient is 

C(F - A) = -CA1 (Same as Eq. G-41) 

9. The thrust minus drag coefficient is 

C(F - D) = -CDS1 (Same as Eq. G-42) 
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G-15 


id. The ratio of thrust minus axial force to ideal thrust is 

(F - A)/FI = C(F - A)/CFI (Same as Eq. G-43) 

11. The ratio of thrust minus drag to ideal thrust is 

(F - D)/FI = C(F - D)/CFI (Same as Eq. G-44) 

12. The thrust minus nozzle drag coefficient is 

C(F - DNOZ) = C(F - D) + (CDF AFT - CDFNOZ) (Eq. G-48) 

13. The ratio of thrust minus nozzle drag to ideal thrust is 

(F - DNOZ)/FI = C(F - DNOZ)/CFI (Eq. G-49) 

14. The ratio of thrust to ideal thrust is 

F/FI = (C(F - DNOZ) + CDFNOZ + CDBUn) /CFI (Eq. G-50) 

G. Two-Balance/Afterbody Metric IFAF1 = 1 
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1. The jet axial force coefficient is 
CFJ = CA2 - CA1 

2. The jet normal force coefficient is 
CNJ = CN1 - CN2 

3. The jet pitching moment coefficient is 
CMJ = CMY1 - CMY2 

4. The resultant thrust coefficient about the pitch axis is 
CT = VcFJ 2 + CNJ 2 

5. The thrust vector angle about the pitch axis is 
DELTA = TAN' 1 (CNJ/CFJ) 

6. The thrust removed lift coefficient is 
CLAERO = CLS2 

7. The thrust removed drag coefficient is 
CDAERO = CDS2 

8. The thrust removed pitching moment coefficient is 
CMAERO = CYMS2 

9. The thrust minus axial force coefficient is 
C(F - A) = -CA1 


(Eq. G-51) 


(Eq. G-52) 


(Eq. G-53) 


(Eq. G-54) 


(Eq. G-55) 


(Eq. G-56) 


(Eq. G-57) 


(Eq. G-58) 


(Same as Eq. G-41) 
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10. The thrust minus drag coefficient is 

C(F - D) = -CDS1 (Same as Eq. G-42) 

11. The ratio of thrust minus axial force to ideal thrust is 

(F - A)/FI = C(F - A)/CFI (Same as Eq. G-43) 

12. The ratio of thrust minus drag to ideal thrust is 

(F - D)/FI = C(F - D)/CFI (Same as Eq. G-44) 

13. The thrust minus nozzle drag coefficient is 

C(F - DNOZ) = C(F - D) + (CDFAFT - CDFNOZ) + CDSUN (Eq. G-59) 

14. The ratio of thrust minus nozzle drag to ideal thrust is 

(F - DNOZ)/FI = C(F - DNOZ)/CFI (Same as Eq. G-49) 

15. The ratio of thrust to ideal thrust is 

F/FI = CFJ/CFI (Eq. G-60) 

16. The jet reaction lift coefficient is 

CLJET = CT (siN(ALPHA + DELTaJ (Eq. G-61) 

17. The jet reaction axial force coefficient is 

CFJET = CT (C0S(ALPHA + DELTA)) (Eq. G-62) 
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H. Two-Balance/Afterbody Metric IFAF2 - 1 



NOTE: Wings or calls usually attached 
to balance 2 


1. The thrust removed lift coefficient is 

CLAERO = CLS2 (Same as Eq. G-56) 

2. The thrust removed drag coefficient is 

CDAERO = CDS 2 (Same as Eq. G-57) 

3. The thrust removed pitching moment coefficient is 

CMAERO = CYMS2 (Same as Eq. G-58) 

4. The thrust minus axial force coefficient is 

C(F - A) = -CA1 (Same as Eq. G-41) 

5. The thrust minus drag coefficient is 


C(F - D) = -CDS1 


(Same as Eq. G-42) 
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6. The ratio of thrust minus axial force to ideal thrust 
(F - A)/FI = C(F - A)/CFI 

7. The ratio of thrust minus drag to ideal thrust is 
(F - D)/FI = C(F - D)/CFI 

8. The thrust minus nozzle drag coefficient is 
C(F - DNOZ) = C(F - D) + CDS2 

9. The jet axial force coefficient is 

CFJ = C(F - DNOZ) + (CDFNOZ + CDBUN) 

10. The jet normal force coefficient is 
CNJ = CN1 - CN2 - CNBUN 

11. The jet pitching moment coefficient is 
CMJ = CMY1 - CMY2 - CPMBUN 

12. The resultant thrust coefficient about the pitch axis 
CT = VcFJ 2 + CNJ 2 

13. The thrust vector angle about the pitch axis is 
DELTA = TAN -1 (CNJ/CFJ) 

14. The jet reaction lift coefficient is 
CLJET = CT (SIN(ALPHA + DELTA)) 


is 

(Same as Eq. G-43) 

(Same as Eq. G-44) 

(Eq. G-63) 

(Eq. G-64) 

(Eq. G-65) 

(Eq. G-66) 
is 

(Same as Eq. G-54) 
(Same as Eq. G-55) 
(Same as Eq. G-61) 
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15. The jet reaction axial force coefficient is 


CFJET = CT (C0S(ALPHA + DELTA)) (Same as Eq. G-62) 

Two-Balance/Afterbody Metric IFAFN1 = 1 



1. The jet axial force coefficient is 

CFJ = -CA2 (Eq. G-67) 

2. The jet normal force coefficient is 

CNJ = CN2 (Eq. G-68) 

3. The jet pitching moment coefficient is 

CMJ = CMY2 (Eq. G-69) 

4. The resultant thrust coefficient about the pitch axis is 
CT = \/cFJ 1 2 3 4 + CNJ 2 


(Same as Eq. G-54) 
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5. The thrust vector angle about the pitch axis is 
DELTA = TAN -1 (CNJ/CFJ) 

6. The jet reaction lift coefficient is 
CLJET = CT (SIN(ALPHA + DELTA)] 

7. The jet reaction axial force coefficient is 
CFJET = CT (C0S(ALPHA + DELTA)) 

8. The thrust minus axial force coefficient is 
C(F - A) = -CA1 

9. The thrust minus drag coefficient is 
C(F - D) = -CDS1 

10. The ratio of thrust minus axial force to ideal thrust 
(F - A)/FI = C(F - A)/CFI 

11. The ratio of thrust minus drag to ideal thrust is 
(F - D)/FI = C(F - D)/CFI 

12. The thrust removed lift coefficient is 
CLAERO = CLS1 - CLS2 

13. The thrust removed drag coefficient is 
CDAERO = CDS1 - CDS2 


(Same as Eq. G-55) 


(Same as Eq. G-61) 


(Same as Eq. G-62) 


(Same as Eq. G-41) 


(Same as Eq. G-42 


(Same as Eq. G-43) 


(Same as Eq. G-44) 


(Eq. G-70) 


(Eq. G-71) 
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G-22 


14. The thrust removed pitching moment coefficient is 

CMAERO = CMYS1 - CMYS2 (Eq. G-72) 

15. The thrust minus nozzle drag coefficient is 

C(F - DNOZ) = C(F - D) + (CDF AFT - CDFNOZ) + CDSUN 

(Same as Eq. G-59) 


J. Two-Balance/Afterbody Metric IFAFN2 = 1 



1. The jet axial force coefficient is 


CFJ = (CDFNOZ + CFBUN) - CA2 
2. The jet normal force coefficient is 


CNJ = CN2 - CNBUN 


(Eq. G-73) 


(Eq. G-74) 


3. The jet pitching moment coefficient is 


CMJ = CMY2 - CP MB UN 


(Eq. G-75) 


4. 


The resultant thrust coefficient about the 
CT = /cFJ 2 + CNJ 2 


pitch axis is 

(Same as Eq. G-54) 
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G-23 

5. The thrust vector about the pitch axis is 

DELTA = TAN -1 (CNJ/CFJ) (Same as Eq. G-55) 

6. The jet reaction lift coefficient is 

CLJET = CT * [SIN (ALPHA + DELTA)) (Same as Eq. G-61) 

7. The jet reaction axial force coefficient is 

CFJET = CT * (C0S(ALPHA + DELTA)) (Same as Eq. G-62) 

8. The thrust minus axial force coefficient is 

C(F - A) = -CA1 (Same as Eq. G-41) 

9. The thrust minus drag coefficient is 

C(F - D) = -CDS1 (Same as Eq. G-42) 

10. The ratio of thrust minus axial force to ideal thrust is 

(F - A)/FI = C(F - A)/CFI (Same as Eq. G-43) 

11. The ratio of thrust minus drag to ideal thrust is 

(F - D)/FI = C(F - D)/CFI (Same as Eq. G-44) 

12. The thrust removed lift coefficient is 

CLAERO = CLS1 - CLS2 (Same as Eq. G-70) 

13. The thrust removed drag coefficient is 

CDAERO = CDS1 - CDS2 (Same as Eq. G-71) 
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14. The thrust removed pitching moment coefficient is 


CMAERO = CMYS1 - CMYS2 
15. The thrust minus nozzle drag coefficient is 
C(F-DNOZ) = -CDS2 

Single Balance. Thrust Removal All Components IFAF 

1. The thrust removed normal force coefficient is 

CNAERO = CN1 - CNJC 

2. The thrust removed axial force coefficient is 

CAAERO = CA1 + CFJC 

3. The thrust removed pitching moment coefficient is 

CMAERO = CM Y1 - CMJC 

4. The thrust removed rolling moment coefficient is 

CRAERO = CMX1 - CRJC 

5. The thrust removed yawing moment coefficient is 

CYAERO = CMZ1 - CYJC 

6. The thrust removed side force coefficient is 


(Same as Eq. G-72) 


(Eq. G-76) 


(Eq. G-77) 


(Eq. G-78) 


(Eq. G-79) 


(Eq. G-80) 


(Eq. G-81) 


CSAERO = CY1 - CSJC 


(Eq. G-82) 


7. The thrust removed lift coefficient is 


G-25 


CLAERO = CNAERO * COS(ALPHA) - CAAERO * SIN(ALPHA) 

(Eq. G-83) 


8. The thrust removed drag coefficient is 

CDAERO = CAAERO * COS(ALPHA) + CNAERO * SIN(ALPHA) 

(Eq. G-84) 

9. The resultant thrust coefficient about the pitch axis is 

CT = ^/cNJC^ + CFJC^ (Same as Eq. G-34) 


10. The thrust vector angle about the pitch axis is 

DELTA = TAN' 1 (CNJC/CFJC) (Same as Eq. G-35) 

11. The jet reaction lift coefficient is 


CLJET = CNJC * COS (ALPHA) + CFJC * SIN (ALPHA) 

12. The jet reaction axial force coefficient is 

CFJET = CFJC * COS (ALPHA) - CNJC * SIN (ALPHA) 

13. The jet reaction side force coefficient is 

CSJET = CSJC 


14. The jet reaction pitching moment coefficient is 
CM JET + CMJC 
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15. The jet reaction rolling moment coefficient is 

CRJET = CRJC * COS (ALPHA) + CYJC * SIN (ALPHA) 

16. The jet reaction yawing moment coefficient is 

CYJET = CYJC * COS (ALPHA) - CRJC * SIN (ALPHA) 

17. The splay angle is 

SPLAY1 = TAN" 1 (CSJC/CNJC) 

18. The thrust minus axial force coefficient is 

C(F - A) = -CA1 (Same as Eq. G-41) 

19. The thrust minus drag coefficient is 

C(F - D) = -CDS1 (Same as Eq. G-42) 

20. The ratio of thrust minus axial force to ideal thrust is 

(F - A) /FI = C(F - A)/CFI (Eq. G-85) 

21. The ratio of thrust minus drag to the ideal thrust is 

(F - D)/FI = C(F - D)/CFI (Eq. G-86) 

22. The ratio of thrust to the ideal thrust is 


F/FI = CFJC/CFI 


(Same as Eq. G-60) 


G-27 


23. The resultant thrust coefficient about the yaw axis is 

CTY = n/cSJC 2 + CFJC 2 

24. The thrust vector angle about the yaw axis is 

DELTAY = TAN -1 (CSJC/CFJC) 

25. The total resultant thrust coefficient is 

CTT = \/CNJC 2 + CFJC 2 + CSJC 2 


(Eq. G-87) 


(Eq. G-88) 


(Eq. G-89) 


26. The total resultant thrust ratio is 


FT/FI = CTT/CFI 


(Eq. G-90) 


27. The thrust minus nozzle axial force coefficient is 

C(F-ANOZ) = C(F-A) + (CDFAFT - CDFNOZ) + CFSUN 


28. The ratio of thrust minus nozzle axial force to the ideal thrust is 
F-ANOZ/FI = C(F-ANOZ)/CFI 


29. The thrust minus nozzle drag coefficient is 

C(F-DNOZ) = C(F-D) + (CDFAFT - CDFNOZ) + CDSUN 

(Same as Eq. G-59) 

30. The ratio of thrust minus nozzle drag to the ideal thrust is 

(F-DNOZ)/FI = C(F-DNOZ)/CFI 
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31. The thrust coefficient (from balance and pressures) is 


CF = (F-ANOZ) + CDFNOZ + CFBUN 

32. The ratio of thrust (from balance and pressures) to the ideal thrust is 
CF/CFI = CF/CFI 
When 1FAFN = 1 

1. The thrust removed stability axis lift coefficient is 

CLSAER = CLAERO 

2. The thrust removed stability axis drag coefficient is 

CDSAER = CDAERO 

3. The thrust removed stability axis side force coefficient is 

CSSAER = CSAERO 

4. The thrust removed stability axis pitching moment coefficient is 

CMSAER = CMAERO 

5. The thrust removed stability axis rolling moment coefficient is 

CRSAER = CRAERO COS ALPHA + CYAERO SIN ALPHA 

6. The thrust removed stability axis yawing moment coefficient is 

CYSAER = CYAERO COS ALPHA - CRAERO SIN ALPHA 



7. The thrust removed wind axis drag coefficient is 

CDWAER = CDSAER COS BETA - CSSAER SIN BETA 

8. The thrust removed wind axis side force coefficient is 

CSWAER = CSSAER COS BETA + CDSAER SIN BETA 

9. The thrust removed wind axis lift coefficient is 

CLWAER = CLSAER 

10. The thrust removed wind axis rolling moment coefficient is 

CRWAER = CRSAER COS BETA + CMSAER SIN BETA 

11. The thrust removed wind axis pitching moment coefficient is 

CMWAER = CMSAER COS BETA - CRSAER SIN BETA 

12. The thrust removed wind axis yawing moment coefficient is 

CYWAER = CYSAER 

Bifurcate Support Mode Two Balance/Afterbody Metric IDN = 1 
1. The axial force coefficient is modified 
CA1 = CA1 - 0.0004 


2. The drag coefficient in the stability axis is modified 


CDS1 = CDS1 - 0.0004 


The thrust removed normal force coefficient is 


CNAERO = CN1 - CNJC (Same as Eq. G-77) 

The thrust removed axial force coefficient is 

CAAERO = CA1 + CFJC (Same as Eq. G-78) 

The thrust removed pitching moment coefficient is 

CMAERO = CMY1 - CMJC (Same as Eq. G-79) 

The thrust removed rolling moment coefficient is 

CRAERO = CMX1 - CRJC (Same as Eq. G-80) 

The thrust removed yawing moment coefficient is 

CYAERO = CMZ1 - CYJC (Same as Eq. G-81) 

The thrust removed side force coefficient is 

CSAERO = CY1 - CSJC (Same as Eq. G-82) 

The thrust removed lift coefficient is 

CLAERO = CNAERO * COS (ALPHA) - CAAERO * SIN (ALPHA) 

(Same as Eq. G-83) 

The thrust removed drag coefficient is 

CDAERO = CAAERO * COS (ALPHA) + CNAERO * SIN (ALPHA) 

(Same as Eq. G-84) 


G-31 


11. The computed resultant thrust about the pitch axis is 

CT = v/cNJC 2 + CFJC 2 (Same as Eq. G-34) 

12. The computed thrust vector angle about the pitch axis is 

DELTA = TAN -1 (CNJC/CFJC) (Same as Eq. G-35) 


13. The jet reaction lift coefficient is 

CLJET = CtJsIN (ALPHA + DELTA)J (Same as Eq. G-36) 

14. The jet reaction axial force coefficient is 

CFJET = Ct|cOS (ALPHA + DELTA)] (Same as Eq. G-37) 

15. The thrust minus axial force coefficient is 

C(F-A) = -CA1 (Same as Eq. G-41) 

16. The thrust minus drag coefficent is 


C(F-D) = -CDS1 


(Same as Eq. G-42) 


17. The ratio of thrust minus axial force to ideal thrust is 


(F-A)/FI = C(F - AF)/CFI (Same as Eq. G-43) 

18. The ratio of thrust minus drag to ideal thrust is 

(F - D)/FI = C(F - D)/CFI (Same as Eq. G-44) 
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The ratio of thrust to ideal thrust is 


F/FI = CFJC/CFI (Same as Eq. G-45) 

The computed resultant thrust coefficient about the yaw axis is 

CTY = JcYJC 2 + CFJC 2 (Same as Eq. G-87) 

The computed thrust vector angle about the yaw axis is 

DELTAY = TAN -1 (CYJC/CFJC) (Same as Eq. G-88) 

The computed resultant thrust coefficient is 

CTT = n/cNJC 2 + CFJC 2 + CYJC 2 (Same as Eq. G-89) 

The total resultant thrust ratio is 

FT/FI = CTT/CFI (Same as Eq. G-90) 

The thrust minus nozzle drag coefficient is 
C(D-FNOZ)/FI = CDS2 - CDS1 
The ratio of thrust minus nozzle drag to ideal thrust is 
(F - FNOZ)/FI = C(D - FNOZ)/CFI 
The nozzle lift plus thrust coefficient is 


CLNOZT = CLS1 - CLS2 
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27. The nozzle pitching moment plus lift coefficient is 

CMNOZT = CMYS1 - CMYS2 

28. The nozzle drag coefficient is 

CDNOZ = CDAERO - CDS2 

29. The thrust removed nozzle lift coefficient is 

CLNOZ = CLAERO - CLS2 

30. The thrust removed nozzle pitching moment coefficient is 

CMNOZ = CMAERO - CMYS2 
N. Other Options 

1. ID - Engineer's option 

a. If ID = 1, the engineer may write his own option with the following 
restrictions: 

(1) Names must be identical to those already used. 

(2) No more terms may be added to the output. 
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MODULE H TURBOPROP OPTIONS 


H-l 


SYMBOL 

NOMENCLATURE 

AD(J) 

Area at rake in exhaust duct J, sq. in. 

AE(J) 

Exhaust area for exhaust duct J, sq. in. 

ALPHAP(J) 

Angle of attack at propeller J, degrees. 

ARATIO(J) 

Area ratio for motor J. 

AT(J) 

Throat area of exhaust duct J, sq. in. 

CDP(J) 

Propeller drag coefficient for motor J. 

CDTP 

Total propeller drag coefficient. 

CE 

Exhaust sonic velocity, feet per second. 

CO 

Sonic velocity, feet per second. 

CHPROP(J) 

Chord length at 75% radius of propeller J. 

CMPROP(J) 

Pitching moment coefficient of propeller J. 

CNPROP(J) 

Normal force coefficient of propeller J. 

CPPROP(J) 

Power coefficient of propeller J. 

CTPROP(J) 

Thrust coefficient of propeller J. 

DIAP(J) 

Diameter of propeller J, feet. 

ETA(J) 

Efficiency for motor J, per cent. 

ETAP(J) 

Efficiency for propeller J, per cent. 

FAPT 

Total system thrust in streamwise direction, lbs. 

FTE(J) 

Propeller thrust plus jet thrust due to exhaust flow for 
motor J, lbs. 

FTGE(J) 

Total system thrust of motor J, lbs. 

JBINC 

Increment on engine number to match prop balance number. 

JP(J) 

Advance ratio of propeller J. 

KBINC 

Indicates that balance 3 deck contains more than one set of 
balance constants. (KBINC = 1). 

KPINM(I,J) 

Constant for input drive pressure tap I and motor J, (must 
be 0.0 or 1.0). 

KPOUTM(I,J) 

Constant for output drive pressure tap I and motor J, (must 
be 0.0 or 1.0). 

KPST(I,J) 

Constant for rake static pressure tap I and motor J, (must 
be 0.0 or 1.0). 

KPW 

Power coefficient constant. 
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H-2 

SYMBOL 

KTINM(I,J) 

KTOUTM(I,J) 

MD(J) 

ME(J) 

MTIP(J) 

NPROP(J) 

NSAME(J) 

PDRIVE 

PE(J) 

PHIANG 

PINM(I,J) 

PITCH(J) 

POUTM(I,J) 

PR/PTR 

PST(I,J) 

PSTATC(J) 

PW1(J) 

PW2(J) 

RHO 

RPS 

TDRIVE(J) 

TE(J) 

TINM(I,J) 

TOUTM(I,J) 

TSPROP(J) 

TTO 

VE(J) 


NOMENCLATURE 

Constant for input motor temperature tap I and motor J, 

(must be 0.0 or 1.0). 

Constant for output motor temperature tap I and motor J, 

(must be 0.0 or 1.0). • 

Rake mach number for motor J. 

■ 

Exhaust mach number for motor J. 

Mach number of propeller tip J. 

Revolutions per second of propeller J. 

Constant of propeller J set equal to 0.0 or 1.0 
Pressure drop through air turbine motor, lbs/sq. in. 

Exhaust static pressure for motor J, lbs/sq. in. 

Angle between forward and rotational velocities, degrees. 

Motor input static pressure for motor J and pressure tap I, 
lbs/sq. in. 

Measured value of geometric pitch of propeller J. 

Motor output static pressure for motor J and pressure tap I, 
lbs/sq. in. 

Ratio of static to total pressures. 

Static pressure for motor J and pressure tap I at rake, 
lbs/sq. in. 

Average static pressure at rake for motor J, lbs/sq. in. 

Horsepower output by motor J with ideal gas calculations, 

HP. 

Horsepower calculated using Isentropic equation multiplied 
by efficiency for motor J, HP. 

Density of free-stream air, slugs per cubic feet. 

Engine's revolutions per second. 

Temperature differential across the air turbine motor J, F. 

Exhaust temperature for motor J, F. 

Input temperature for motor J and temperature tap I, F. 

Output temperature for motor J and temperature tap I, °F. 

Rotational tip speed of propeller J, feet per second. 

Tunnel static temperature, °F. 

Exhaust velocity in motor J, feet per second. 
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H-3 


SYMBOL 

VO 

VRES(J) 

VRN 


NOMENCLATURE 

Free-stream velocity, feet per second. 

Total velocity of propeller tip J, feet per second. 

Total velocity at 75% of propeller radius (for Reynolds 
number), feet per second. 
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Module H 

Turboprop Options 


Introduction 

1. Module B with its constants must be run first. All constants are to be 
initialized to a value of zero. The project engineer must supply only those 
constants which are required for those quantities to be computed. In addition, 
by logical use of combinations of these constants, several options are available 
to the project engineer. 

2. Set NSAME(J) = 1 if POUTM(I,J) = PST(I,J), and TOUTM(I,J) = TTJ(I,J) 

where J = engine number 
I = probe number 

3. Set the constant, KTINM(I,J), equal to 1.0 for the temperature measuring probe. 
If the temperature probe is defective or does not exist, set the constant equal 
to 0.0. Use only a maximum of six probes per engine. 

4. The meaning of the values of KTOUTM(I,J) is the same as KTINM(I,J). 

5. Set the constant, KPINM(I,J) equal to 1.0 for the pressure measuring probe. If 

the probe is defective or does not exist, set the constant equal to 0.0. The 

pressure probes may be weighted using this constants, if desired. Use only a 
maximum of 12 probes per engine. 

6. The meaning of the values of KPOUTM(I,J) is the same as KPINM(I,J). 

7. Set the constant, KPST(I,J), equal to 1.0 for the pressure measuring probe. If 

the probe is defective or does not exist, set the constant equal to 0.0. The 

pressure probes may be weighted using this constant, if desired. Use only a 
maximum of 12 probes per engine. 
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8. AE(J) is equal to AT(J) for a converging nozzle. Both constants are required. 
Values of AT(J) come from Module B. 


B. Test for Air Turbine Simulator 

1. The constant required from the project engineer input at Module B is NUMENG 
(0 to 4). 

If NUMENG = 0, skip module H. 

C. Compute Common Constant 


1. The constants required from the project engineer input at Module B are GAMJ 
and RJ. 


KJ1 = 


( 


2 _) 

GAMJ + 1 / 


GAMJ + 1 
2(GAMJ - 1) 


yGAMJ 


x 32.174 
RJ 


(Same as Eq. B-l) 


GAMJ x 64.348 
(GAMJ - 1)RJ 


(Same as Eq. B-2) 


KJ3 = 


A 


2(GAMJ) (RJ) 


(GAMJ - 1)32.174 


(Same as Eq. B-3) 


KJ4 = 


GAMJ - 1 
GAMJ 


(Same as Eq. B-4) 


KJ5 = 


1 _ 

GAMJ 


(Same as Eq. B-5) 
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To continue, the equations are given to show calculations for other constants. 


a. The static temperature is 

TO = (TTO + 459.67)/(1.0 + 0.2 * MACH 2 ) 

b. The free-stream density is 

RHO = PO * 144.0/(1716.4829 * TO) 

c. The viscosity is 

XMU = (2.270 * 10" 8 * TO * VTO/(TO + 198.6) 

d. The free-stream velocity of sound is 
CO = 49.021179 *VTO 

e. The free-stream velocity is 
VO = CO * MACH 

Individual Engine Measurements 

1. This module provides the computations for four separate 
following instrumentation in each engine. 

a. Input pressure to engine 

*b. Output pressure of engine 

c. Input temperature to engine 

*d. Output temperature of engine 


(Eq. H-l) 


(Eq. H-2) 


(Eq. H-3) 


(Eq. H-4) 


(Eq. H-5) 


engines with the 


♦ May be replaced with rake measurements. 



e. Static exhaust pressure at rake 

f. Revolutions per second indicator 

g. Geometric pitch of propeller 
Propeller Coefficient Calculations 

1. The tip speed of propeller is 

TSPROP(J) = 3.14159 * DIAP(J) * NPROP(J) (Eq. H-6) 

2. The advance ratio of propeller is 

JP(J) = VO/(NPROP(J) * DIAP(J)) (Eq. H-7) 

3. The angle of attack of the propeller is the geometric pitch of the propeller at 
the 3/4 chord, in degrees, minus the resultant angle between the free-stream 
velocity and rotational velocity. 

ALPHAP(J) = PITCH(J) - PHIANG (Eq. H-8) 

where 

PHIANG = TAN" 1 (VO/VROT) (Eq. H-9) 

and 


VROT = 3/4 TSPROP(J) 


(Eq. H-10) 
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The Reynolds number for the propeller is calculated at the 3/4 chord. 
RNPROP(J) = VRN * RHO * CHPROP(J)/XMU (Eq. H-ll) 

where 

VRN = Resultant velocity at the 3/4 chord 
= V (VROT 2 + VO 2 ) 

XMU = Free-stream air viscosity calculated by Ames table equation, based on 
tunnel air static temperature 

The Mach number of the propeller tip is 

MTIP(J) = VRES/CO (Eq. H-12) 


where 


VRES = 



+ TSPROP(J) 2 


Calculate the thrust coefficient of the propeller and hub using 
SCALE = RHO * NPROP(J) 2 * DIAP(J) 4 
CTPROP(J) = FAREF1/SCALE 


(Eq. H-13) 

(Eq. H-14) 
(Eq. H-15) 


where 

FAREF1 comes from Equation D-75. 

Calculate the normal force coefficient of the propeller and hub using 


CNPROP(J) = FNREF1/SCALE 


(Eq. H-16) 



H-10 


where 


FNREF1 comes from Equation D-75. 

8. Calculated the pitching moment coefficient of the propeller and hub using 

CMPROP(J) = MYREF1/(SCALE * DIAP(J) * 12.0) (Eq. H-17) 

where 

MYREF1 comes from Equation D-75. 

9. If NSAME(J) equals 1, then 

POUTM(I,J) = PST(I,J) (Eq. H-18) 

and 


TOUTM(I,J) = TTJ(I,J) 


(Eq. H-19) 


10. Calculations for the power coefficient of the propeller and hub are: 


a. Turbine inlet temperature 

£lTNM(I,J) * KTINM(I,J) 


TIN(J) = 


£KTINM(I,J) 


(Eq. H-20) 


b. Turbine outlet temperature 

?TOUTM(I,J) * KTOUTM(I,J) 
TOUT(J) = — 

£ktoutm(i,j) 


(Eq. H-21) 


c. Turbine inlet pressure 

PIN(J) = yPINM(I,J) * KPINM(I,J) 
£kPINM(I,J) 


(Eq. H-22) 
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H-ll 


d. 


Turbine outlet pressure 

P0UT(J) = ;POUTM(l,J) ♦ KPOUTM(l,J) 

l KPOUTM(I,J) 


(Eq. H-23) 


e. The drive pressure across the air turbine engine is 


PDRIVE(J) = PIN(J) - POUT(J) 


(Eq. H-24) 


f. The drive temperature across the air turbine engine is 

TDRIVE(J) = TIN(J) - TOUT(J) (Eq. H-25) 


g. The engine's revolutions per second are 


RPS = NPROP(J)/ V(TIN(J) + 459.67)7518.7 (Eq. H-26) 


h. Calculate the horsepower output from the air turbine engine using 


PVV1(J) = (6006.0 * (WPENG(J)/32.174) * TDRIVE(J)/550) * (Eq. H-27) 

((KPW13 * PIN(J) + KPW12 * RPS + KPW11) * RPS + KPW10) 

PW2(J) = (6006.0 * (WPENG(J)/32.174) 

* (TIN(J) + 459.67) 

* (1.0 - (POUT(J)/PIN(J)) 2/7 ) 

* ETA(J)/550 (Eq. H-28) 

where 


ETA(J) is determined by linear interpolation from a table. 
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i. The power coefficient of the propeller and hub is 

CPPROP(J) = PW/(RHO * NPROP(J) 3 * DIAP(J) 5 ) (Eq. H-29) 

PW = PW2(J) if KPW = 0 
PW = PW1(J) if KPW = 1 

j. The propeller efficiency is 

ETAP(J) = CTPROP(J) * JP(J)/CPPROP(J) (Eq. H-30) 


Exhaust Calculations 


1. Calculate exhaust duct Mach number (rake position) using 


a. 


Duct static pressure 


PSTATIC(J) = 


l PST(I, J) * KPST(I,J) 
l KPST(I,J) 


(Eq. H-31) 


b. The pressure ratio at the duct rake is 

PR/PTR(J) = PSTATIC(J)/PTEN G(J) (Eq. H-32) 

c. If PR/PTR(J) = PSTATIC(J) < .5283, use the Newton Raphson method for 
MD(J). 

«*«« Vt> ( 7 * T )2r *) (UkOt) ' W ** 


d. If PR/PTR(J)> .5283, use this calculation of subsonic duct Mach numbers 
for MD(J). 

MD(J) = (PR/PTR(J)) ~ 2/7 - 5 


(Eq. H-34) 



The ratio of A* to Area at the rake position and at the exit is 


H-13 


a. 


Calculate A* /A at the rake station using 


ASTR/A = (1.728 * MD(J)) * 



MD(J) 2 ^ 


b. Calculate A* /A of the exhaust exit using 

ARATIO(J) = ASTR/A * AD(J)/AE(J) 


Calculate the exhaust Mach number at the exit using an 
the formula 


ME(J) = 


* (ARATIO(J)) * 



ME(J) 

5 



The exhaust static temperature calculation is 


TE = (TTENG(J) + 459.67) * 



i ME(J) 2 



where 


TTENG(J) comes from Equation B-9. 


The exhaust sonic velocity is 
CE = 49.021179 * VtE 


The exhaust velocity is 
VE(J) = ME(J) * CE 


The exhaust static pressure is 


(Eq. H-35) 
(Eq. H-36) 

iteration technique on 
(Eq. H-37) 

(Eq. H-38) 


(Eq. H-39) 
(Eq. H-40) 
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H-14 


PE(J) = PTENG(J) * (l + ME( 5 J) ) 


(Eq. H-41) 


-7/2 


8. The total propeller pitching moment is 


MYPT = V MY. 

L 1 


9. The total propeller normal force is 

NUMENG 

FNPT = l NF. 

i = 1 1 


10. The total propeller axial force is 


FAPT = I AF. 

i = l 1 


(Eq. H-42) 


(Eq. H-43) 


(Eq. H-44) 


11. The axial force coefficient in the body axis with propeller and jet thrust 
removed is 

CAPRS = CAAERO + FAPT (Eq. H-45) 

12. The drag coefficient in the stability axis with propeller and jet thrust removed 
is 

CDPRS = CAPRS COS ALPHA + CNAERO SIN ALPHA (Eq. H-46) 

13. The side force coefficient in the stability axis with propeller and jet thrust 
removed is 

CSPRS = CSAERO (Eq. H-47) 
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The lift coefficient in the stability axis with propeller and jet thrust removed is 


CLPRS = CNAERO COS ALPHA - CAPRS SIN ALPHA (Eq. H-48) 

The rolling moment coefficient in the stability axis with propeller and jet thrust 
removed is 

CRPRS = CRAERO COS ALPHA + CYAERO SIN ALPHA (Eq. H-49) 

The pitching moment coefficient in the stability axis with propeller and jet 
exhaust removed is 

CMPRS = CMAERO (Eq. H-50) 

The yawing moment coefficient in the stability axis with propeller and jet 
exhaust thrust removed is 

CYPRS = CYAERO COS ALPHA - CRAERO SIN ALPHA (Eq. H-51) 

The drag coefficient in the wind axis with propeller and jet thrust removed is 

CDPRW = CDPRS COS BETA - CSPRS SIN BETA (Eq. H-52) 

The side force coefficient in the wind axis with propeller and jet thrust 
removed is 

CDPRW = CSPRS COS BETA + CDPRS SIN BETA (Eq. H-53) 

The lift coefficient in the wind axis with propeller and jet exhaust thrust 
removed is 


CLPRW = CLPRS 


(Eq. H-54) 


H-16 


21. The polling moment coefficient in the wind axis with propeller and jet exhaust 
thrust removed is 

CRPRW = CRPRS COS BETA + CMPRS SIN BETA (Eq. H-55) 

22. The pitching moment coefficient in the wind axis with propeller and jet exhaust 
thrust removed is 

CMPRW = CMPRS COS BETA - CRPRS SIN BETA (Eq. H-56) 

23. The yawing moment coefficient in the wind axis with propeller and jet exhaust 
thrust removed is 


CYPRW = CYPRS 


(Eq. H-57) 
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